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DISTRIBUTION  STATEMENT 


ABSTRACT 


An  experimental  and  theoretical  study  was  made  of  a mono- 
pole antenna  mounted  on  a finite  ground  plane  located  above  an  infinite 
ground.  Circular  flat  discs  and  hemispheres  were  used  for  the  finite 
ground  planes.  Experiments  were  conducted  over  a 6 to  1 frequency 
range,  with  the  length  of  monopole  antenna  fixed  at  a quarter  wave- 
length at  the  upper  end  of  the  frequency  band.  The  radii  of  finite 
ground  planes  used  for  experiments  were  generally  less  than  a quarter 
wavelength. 

Input  impedances  were  measured  as  a function  of  frequency 
at  the  base  of  the  antenna,  using  as  variables  the  radius  and  the  locations 
of  the  ground  plane  with  respect  to  an  infinite  ground  below.  Scale 
models  were  used  to  obtain  measurements  of  radiation  patterns  and 
antenna  current  distributions.  Results  of  these  measurements  are 
presented  graphically. 

A theoretical  analysis  was  also  made  of  a monopole  with  a 
hemispherical  ground  plane  on  an  infinite  ground.  Far-zone  electro- 
magnetic fields  were  calculated  as  a function  of  both  the  ground  plane 
radius  and  the  frequency.  Radiation  resistances  were  also  calculated. 

The  results  of  the  study  indicate  that  the  radiation  resistance 
of  an  electrically  short  antenna  may  be  increased  significantly  by 
locating  it  on  a small  ground  plane  above  the  infinite  ground,  rather 
than  directly  on  the  infinite  ground.  This  conclusion  opens  the  way  for 


iii 


a more  efficient  utilization  of  receiving  and  transmitting  antennas  on 
ground-  hised  vehicles.  The  impedance  characteristics  of  the  antenna 
system  art;  such  as  to  facilitate  its  operation  over  a broad  frequency 
range. 
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FOREWORD 


This  report  was  prepared  by  the  Cooley  Electronics  Labora- 
tory of  The  University  of  Michigan  under  United  States  Army  Electronics 
Command  Contract  No.  DA  28-043-AMC-02246(E),  Project  No.  5A6-79191- 
D902-02-24,  "Improved  Antenna  Techniques  Study.  " 

The  research  under  this  contract  consists  in  part  of  an 
investigation  to  develop  highly  efficient  remotely  tuned  impedance  match- 
ing coupling  networks  for  electrically  short  monopoles. 

The  material  reported  herein  represents  a summary  of  a 
theoretical  and  experimental  study  which  was  made  to  determine  the  input 
impedance  and  radiation  characteristics  of  an  electrically  short  monopole 
over  a small  ground  plane  located  at  various  distances  above  natural 
ground. 
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CHAPTER  I 


INTRODUCTION 

Monopole  and  dipole  antennas  ov^r  an  infinitely  large  con- 
ducting ground  plane  are  two  of  the  oldest  types  of  antennas  being  used. 
Their  performance  characteristics  have  been  thoroughly  studied  in  the 
past.  However,  the  exact  behavior  of  the  antenna  over  a conducting 
ground  plane  of  a finite  size,  and  both  antenna  and  finite  ground  plane 
over  infinite  ground  with  finite  conductivity  has  not  been  studied  here- 
tofore. 

1. 1 Statement  of  the  Problem 

The  purpose  of  this  study  is  to  investigate  the  properties 
of  a vertical  monopole  antenna  and  a finite  ground  plane  located  above 
an  infinite  ground  plane.  Particular  emphasis,  in  this  investigation, 
is  placed  upon  the  study  of  antenna  systems  having  circular  ground  planes 
m which  the  dimensions  of  both  the  ground-plane  diameters  and  the  ground- 
piano  heights  above  infinite  ground  are  smaller  than  the  wavelengths 
of  interest. 

The  investigation  of  the  properties  of  this  antenna  system 
provides  a basis  for  a detailed  comparison  with  the  conventional  antenna 
system  of  a monopole  antenna  over  an  infinite  ground  plane  and  for  de- 
termvntion  of  the  relative  advantages  and  limitations  of  this  system. 

In  iition,  the  resuHs  of  this  study  will  be  applied  to  an  investigation 
of  possible  performancp  improvement  for  a monopole  antenna  mounted 
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on  a ground-based  vehicle,  since  there  is  reason  to  believe  that  the  an- 
tenna system  proposed  for  study  approximates  the  vehicular-mounted 
antenna. 

1.  2 Topics  of  Investigation 

Since  the  purpose  of  this  paper  is  to  report  a study  of  the 
vertical  monopole  with  a finite  ground  plane  in  the  presence  of  an  infi- 
nite ground  plane,  the  following  topics  are  explored: 

(1)  The  experimental  measurement  of  the  input  impedance 
of  the  antenna  was  carried  out  over  the  frequency  range  of  5 MHz  to 
30  MHz  (a  6 to  1 band)  and  as  a funct  ion  both  of  the  ground-plane  size 
and  of  the  ground-plane  location,  relative  to  an  infinite  ground. 

(2)  The  experimental  measurement  of  current  distribution 
on  the  antenna  over  the  frequency  range  described  above  was  carried 
out  with  the  ground  plane  size  and  its  location  as  the  variables. 

(3)  The  radiation  pattern  was  measured  and  its  dependence 
on  the  ground  plane  size  and  location  was  studied. 

(4)  A number  of  theoretical  models  was  examined  to  deter- 
mine a suitable  approximation  to  the  given  system.  Results  oi/ained  from 
the  appropriate  theoretical  model  was  then  ompared  with  the  experimen- 
tal measurements. 

One  of  the  theoretical  models  which  was  studied  is  a mono- 
pole placed  on  a semi- spherical  conductor,  as  shewn  in  Fig.  1. 
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Fig.  1. 1.  Theoretical  models 
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1.  3 Review  of  the  Literature 

Many  studies  of  the  properties  of  vertical  monopole  antennas 
mounted  both  on  infinite  and  on  finite  ground  planes  have  previously  been 
reported  in  the  literature. 

Bardeen  (Ref.  1),  in  1930,  studied  the  diffraction  of  a cir- 
cularly symmetrical  electromagnetic  wave  by  a circular  disc  of  infinite 
conductivity.  The  result  of  his  study  has  been  used  in  determining  the 
power  flow  into  the  earth  below  a vertical  antenna  that  is  grounded  by 
a circular  disc  lying  on  the  surface.  In  1945,  Brown  and  Woodward 
(Ref.  2),  in  an  experimental  study,  measured  the  resistance  and  reac- 
tance of  a cylindrical  antenna  operated  against  ground.  The  ground  im- 
mediately below  the  antenna  had  a buried  metallic  screen  of  a finite 
size  whose  top  surface  was  flush  with  the  rest  of  the  natural  ground. 

A particular  emphasis  was  given  to  a study  of  impedance  behavior  as  a 
function  of  the  antenna  length  to  diameter  ratio.  Terminal  conditions 
such  as  capacitance  of  the  base  of  the  antenna  to  ground  were  considered. 

Later,  Meier  and  Summers  (Ref.  2)  (1949)  performed  an 
experimental  study  of  the  impedance  characteristics  of  vertical  antennas 
mounted  on  finite  ground  planes.  Leitner  and  Spence  (Ref.  3)  (1951) 
confirmed  some  of  these  experimental  results  through  the  theoretical 
study  of  a quarter-wavelength  monopole  on  a finite,  circular-disc  ground 
plane.  They  assumed  a sinusoidal  current  distribution  to  exist  on  the 
antenna.  The  results  of  this  study  showed  a marked  dependence  oi  the 
antenna  radiation  resistance  upon  the  diameter  of  the  disc  employed. 
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This  dependence  was  particularly  pronounced  for  small  ground  planes. 

In  1951,  Storer  (Ref.  4)  used  a variational  method  to  readily  obtain  the 
expression  for  the  dependence  of  the  antenna  impedance  on  the  ground 
plane  diameter.  Tang  (Ref.  5),  in  1962,  further  studied  the  radiation 
pattern  of  this  system  by  using  the  Babinet  principle  and  the  Wiener-Hopf 
technique. 

In  each  of  the  above  papers,  the  monopole  was  considered  to  be 
above  a finite  ground  plane  suspended  in  a free  space  with  no  direct  effect 
from  the  infinite  ground  plane.  Solutions  by  Leitner  and  Spence  (Ref.  3) 
are  mainly  for  the  ground  plane  diameter  smaller  than  the  wavelength 
of  interest.  On  the  other  hand,  Tang  and  Storer* s solutions  can  best 
be  applied  to  the  ground  plane  larger  than  ten  wavelengths  in  diameter. 

Wait  and  Pope  (Ref.  6)  in  1955,  studied  the  input  resistance 
of  low-frequency,  monopole  antennas  with  the  radial-wire  earth  system. 
The  conducting  wires  were  placed  radially  over  the  infinite  natural  ground 
and  the  changes  in  the  input  resistance  AR  from  the  input  resistance  of 
a monopole  with  an  infinite  conducting  ground  plane  was  calculated. 

Brown  (Ref.  7),  in  1937,  and  Abbott  (Ref.  8),  in  1952,  carried 
out  earlier  studies  of  the  radial  wire  system. 

1.  4 Report  Organization 

(1)  Chapter  2 contains  the  development  of  the  experimental 
procedures  for  the  input  impedance  measurement  of  the  combined  antenna 
and  for  the  ground  plane.  The  results  of  both  actual  and  scale  model 
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impedance  measurements  are  included.  Some  discussion  of  the  losses 
involved  in  the  real  part  of  the  impedance  is  included. 

(2)  Chapter  3 contains  a discussion  of  the  technique  for 
making  current  measurements  on  the  antennas.  Results  of  the  measure- 
ments are  shown  in  graphic  form.  Also,  a comparison  is  made  with 

the  sinusoidal  current  distribution  that  is  ordinarily  assumed  for  a mono- 
pole antenna  over  an  infinite  conducting  ground  plane. 

(3)  Chapter  4 contains  a discussion  of  the  technique  of  ra- 
diation-pattern measurements.  The  pertinent  assumptions  made  in  ob- 
taining these  patterns  are  discussed.  The  results  of  the  experiment  are 
displayed  graphically  in  polar  plots  that  show  the  dependence  of  side-lobe 
levels  on  the  ground  plane  size  and  its  location. 

(4)  Chapter  5 contains  an  analysis  of  the  theoretical  model 
outlined  in  the  topics  of  investigation.  The  far-field  electric  and  mag- 
netic fields  are  calculated  for  a given  current  distribution.  The  radiation 
resistances  are  also  calculated  for  different  antennas  and  the  results  are 
compared  with  the  experimental  ones.  The  radiation  patterns  obtained 
from  this  theory  are  also  compared  with  the  experimental  results. 

(5)  Chapter  6 contains  a summary  and  conclusions.  Sug- 
gestions for  further  research  are  presented  and  practical  applications 
usi  ng  the  result  of  this  study  are  discussed. 

The  measurement  procedure  is  described  in  Section  2.  2 
in  which  the  steps  taken  to  ensure  accuracy  in  the  measurement  of  small 
input  resistances  are  given.  The  results  determined  using  this  procedure 
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are  reported  in  some  detail  in  Section  2.  3,  and  compared  with  the  input 
impedance  of  a monopole  on  infinite  ground  plane.  Some  additional 
measurements,  which  are  set  out  in  Section  2. 4,  were  carried  out  with 
scale  models  (1)  of  the  monopole  over  a small  ground  plane  and  both 
above  the  infinite  ground  and  (2)  of  the  monopole  on  a hemispherical 
ground  plane  above  an  infinite  ground.  Experimental  measurements  on 
the  latter  model  were  carried  out  and  were  compared  with  the  theoretical 
results  in  Chapter  5. 


CHAPTER  II 


INPUT  IMPEDANCE  MEASUREMENT 

2. 1 Introduction 

This  chapter  contains  a discussion  of  the  experimental 
work  carried  out  to  measure  the  impedance  of  a monopole  antenna 
that  was  a quarter-wavelength  long  at  30  MHz  and  that  was  located  on 
a finite  ground  that,  in  turn,  was  located  above  natural  ground.  The 
impedance  was  measured  as  a function  both  of  ground  plane  size  and 
of  ground  plane  location. 

Theretofore,  no  extensive  study  of  such  a configuration  has 
been  reported.  The  input  impedance  characteristics  of  an  antenna 
which  is  mounted  on  an  infinite  ground  plane  has  been  reported  by 
Brown  and  Woodward  (Ref.  2)  for  a wide  variety  of  antenna  lengths 
and  diameters.  Both  experimental  and  theoretical  studies  have  been 
conducted  on  the  input  impedance  of  a monopole  antenna  mounted  on  a 
finite  ground  plane,  having  diameters  ranging  down  to  wavelength. 
These  studies  did  not,  however,  consider  the  effect  of  an  infinite 
ground  below  the  finite  ground  plane. 

2.2  Experimental  Measurement  Procedure 

Measurements  were  conducted  with  an  antenna  made  of 
copper  pipe  one-quarter  inch  in  diameter,  which  is  a quarter-wave- 
length long  at  30  MHz.  The  antenna  was  mounted  on  a non-conducting 
variable  height  test  stand,  shown  in  Fig.  2.1.  With  this  stand,  one 
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can  take  measurements  at  heights  of  from  zero  to  5 meters  above 
ground  within  continuous  changes,  and  at  the  ground  plane  diameters 
of  from  zero  to  5 meters. 

In  order  to  observe  the  impedance  variation  with  size  and 
location  of  the  ground  plane  as  measured  in  terms  of  wavelength,  the 
size  and  the  location  of  the  finite  ground  plane  was  varied  in  four 
steps  ranging  from  one- sixteenth  of  a wavelength  to  one- half  wavelength 
at  30  MHz 

The  block  diagram  in  Fig.  2 2 shows  the  test  set-up.  The 
test  equipment  was  placed  just  below  the  ground  plane,  to  minimize 
the  reflection  of  electromagnetic  waves  from  the  equipment.  The 
entire  experiment  was  conducted  in  a flat  field  where  the  nearest 
building  was  at  least  10  wavelengths  away  from  the  test  set-up.  (See 
Fig.  2.  3). 

Four  sizes  of  ground  plane  were  used  in  the  experiments. 

A highly  conducting  5 mil  thick  aluminum  sheet  was  chosen  for  the 
ground  plane  material,  after  considering  the  skin  depth  at  the  frequen- 
cies where  the  measurements  were  taken. 

Because  of  the  provision  for  variable  height  and  diameter, 
the  total  weight  of  the  ground  plane  is  an  important  factor  in  designing 
the  experimental  set-up.  The  depth  of  penetration,  6 , for  aluminum 
at  frequencies  of  5 MHz  and  30  MHz  is  1.  45  mils  and  0.  594  mil,  re- 
spectively. The  thickness  chosen  was  the  smallest, but  was  much 
thicker  than  the  skin  depth  in  aluminum  sheet  ground  plane  A 5 mil 
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Fig.  2.  2.  Block  diagram  showing  the  impedance 
measurement  set-up 


Fig.  2.  4.  A 2.  5 meter  monopole  antenna  on  a 2.  5 
meter  diameter  ground  plane  supported 
by  a styrofoam  sheet 
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thick  roll  of  aluminum  sheeting  is  chosen  for  this  purpose. 

To  support  the  flexible  ground  plane,  a hexagonal  structure 
was  built  and  it  was  covered  with  4 cm  thick  styrofoam  sheet,  which 
has  a dielectric  constant  nearly  equal  to  air.  (Fig.  2. 4).  The  ground 
plane  was  hoisted  up  to  the  proper  height  with  ropes  strung  over  rollers 
located  at  the  top  of  the  six  supporting  posts. 

Initial  measurements  were  carried  out  using  a length  of 
50-ohm  coaxial  line  to  connect  the  base  of  the  antenna  to  the  Impedance 
bridge.  Readings  taken  with  the  bridge  were  referred  to  the  antenna 
base  by  means  of  Smith  charts.  However,  analysis  shows  appreciable 
inaccuracy  to  exist  at  low  frequencies,  where  the  antenna  impedance  ZL  is  high. 
This  inaccuracy  was  due  to  inaccuracies  in  the  impedance  transform- 
ation from  the  bridge  terminals  to  the  antenna  base.  In  transmission 
line  theory,  an  impedance  at  a point  d meters  away  from  the 
load  impedance  ZT  is 


ZL  + j Z0  tan  ~ d 
Z0+iZLtanf  d 


(2.1) 


where  is  the  characteristic  impedance  of  the  line  and  A is  the 
wavelength.  If  the  load  impedance  is  much  greater  thar  the  character- 
istic impedance  Z,  > > Zn  at  any  frequency,  then,  this  expression 


reduces  to 
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Zd  = -jZ0cot-y-d  (2.2) 

which  is  independent  of  the  load  impedance  Z^.  Therefore,  unless  the 
characteristic  impedance  Zq  is  comparable  in  magnitude  to  Z^,  a 
small  change  in  Z^  is  difficult  to  detect  for  large  load  impedances. 
However,  if  d is  chosen  to  be  near  zero  then  the  changes  in  Z^  can 
be  easily  detected  even  if  ZL  > > ZQ. 

From  equation  2. 1 

AZ  . _ Vzo to"2j!rd  AZ 

d 0 (ZQ  + j Zj  tan  d)2  L (2  s) 

and  therefore 


A Z. 


1 Z0  ZT 

cos  2 * + + ~]  sin  26 

1 ZL  Z0 


(2.4) 


where 
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Zird 

\ 


In  general, 


Azd 

zJ 

< 

azl 

ZT 

d 

L 

(2.5) 


Fig.  2.  6.  Location  of  the  impedance  measurement 
bridge  relative  to  the  antenna  base 
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whene/er  Zq  ~ Z^. 

However,  from  the  Eq.  2.  3 when  9 * 0 or  when  d - 0, 


(2.6) 


Equation  2. 6 indicates  that  whenever  the  load  impedance 
ZL  is  much  greater  than  the  characteristic  impedance  of  the  line 
Zq,  the  changes  in  Z L can  be  detected  much  easier  by  choosing  the 
line  length  d = -y-  where  n is  an  integer. 

The  input  impedance  of  a monopole  above  an  infinite  ground 
plane  and  that  resonates  at  a certain  frequency  is  capacitive  at 
frequencies  smaller  than  fQ.  The  magnitude  of  the  capacitive  react- 
ance is  a rapidly  varying  function  of  frequency. 

Therefore,  at  the  lower  end  of  the  frequency  band,  the 
technique  of  using  a line  extension  to  translate  the  impedance  has  a 
large  error;  and,  for  this  reason,  it  was  decided  to  modify  the  ground- 
plane  support  structure  to  accommodate  the  impedance  bridge  directly 
beneath  the  antenna  base.  This  modification  reduced  the  lii  length 
from  the  antenna  to  bridge  to  about  10  cm  and  thereby  required  a 
negligible  impedance  transformation.  The  set-up  is  shown  in  Figs. 

2.  5 and  2.  6. 


Measurements  made  using  this  set-up  are  hampered,  how- 
ever, by  the  limited  range  of  the  bridge.  It  is  capable  of  recording 
only  reactances  of  magnitude  less  than  5, 000 ohms  divided  by  the 
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measurement  frequency  in  megahertz-  The  antenna  capacitive  react- 
ance at  5 MHz  is  greater  than  1000  ohms.  However,  a high-Q  capac- 
itor of  approximately  36  p /if  placed  in  shunt  with  the  antenna  across 
its  base  terminals  reduced  the  overall  reactance  to  the  range  covered 
by  the  bridge,  and  was,  therefore,  used. 

To  obtain  the  input  impedance  of  the  antenna  from  the 
data  obtained  with  the  external  shunt  capacitor,  the  capacitor  was  itself 
carefully  calibrated  over  the  frequency  range  of  5 MHz  to  30  MHz. 
Using  this  equipment,  the  capacitor  impedance,  with  a nominal  value 
of  36  /i  /if,  was  measured  at  every  2.  5 MHz.  The  resistive  component 
and  the  reactive  component  were  carefully  recorded  at  each  frequency. 
These  data  were  later  used  to  obtain  the  actual  input  impedance  of 
the  antenna.  The  resistive  component  of  the  capacitor,  over  the 
frequency  range  of  measurement,  was  always  less  than  one-half  of  an 
ohm.  In  order  to  remove  the  errors  of  manual  calculation,  a computer 
program  was  written  in  which  the  calibration  data  for  the  external 
shunt  caoaciior  and  the  measurement  data  of  the  experiment  were 
fed  in  as  raw  data  and  the  series  resistive  and  reactive  components 
of  the  input  impedance  were  calculated  on  the  IBM  7090  Computer. 

The  calibration  of  the  capacitor  was  done  periodically 
during  the  period  of  this  experiment. 


2.  3 Scale  Model  Impedance  Measurements 

Since  Maxwell's  equations  are  linear  ones,  an  electro- 
magnetic structure  that  has  certain  properties  at  a given  frequency 
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will  have  identical  properties  at  another  frequency  nf,  provided  that 
all  linear  dimensions  of  the  structure  are  scaled  by  the  ratio  1/ n. 
Therefore,  an  antenna  design  which  operates  over  a certain  range  of 
frequencies  can  be  made  to  operate  over  any  other  range  of  frequencies 
without  additional  re-design,  if  an  exact  scaling  of  dimensions  is 
accomplished. 

Aside  from  permitting  one  to  transfer  design  relationships, 
it  is  convenient  to  use  a practical  size  scale  model  for  radiation  pat- 
tern studies.  Scaling  is  a relatively  simple  matter  for  most  types 
of  antennas.  However,  electromagnetic  properties  must  also  be  scaled, 
as  well  as  linear  dimensions.  The  scaling  factors  are  shown  below. 


Length 

1 

1/n 

Frequency 

f 

nf 

Permittivity 

e 

same 

Permeability 

M 

same 

Conductivity 

0 

no 

Table  2. 1.  Scale  factors 

Both  the  length  and  frequency  can  be  scaled  easily,  the 
conductivity  cannot.  However,  the  conductivity  scaling  is  important 
only  through  losses,  and,  since  these  are  small  for  most  antennas, 
the  inability  to  scale  the  conductivity  exactly  is  not  a serious  problem. 

An  exact  scale  model  would  retain  the  exact  radiation 
pattern  and  input  impedance  of  the  full  scale  antenna.  However,  it  is 
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not  always  possible  to  scale  the  model  exactly  t particularly  since 
the  transmission  lines  and  screw  fastenings,  etc. , cannot  be  scaled. 
Slight  discreo.mcies  in  scaling  will  usually  affect  the  impedance 
properties  much  more  than  they  affect  the  radiation  properties.  The 
reason  is  that  any  discrepancies  in  scaling  usually  gives  a rise  in 
near-field  effect  which  are  attenuated  out  at  far-field. 

The  model  used  for  the  antenna  and  ground  plane  studies 
here  is  intended,  therefore,  mainly  for  radiation  pattern  studies.  The 
scale  mode’  designed  to  be  1/40  of  the  original  antenna  system  in 
linear  dimensions.  This  scaling  factor  allows  the  new  antenna  and 
ground  plane  system  to  operate  between  200  MHz  and  1200  MHz, 
replacing  the  original  5 MHz  to  30  MHz  range. 

As  was  pointed  out  previously,  the  principal  difficulty  in 
scaling  lies  with  the  conductivity.  It  is  assumed  here  that  the  copper 
used  for  the  scale  model  ground  plane  and  the  aluminum  ground  plane 
used  for  the  full  size  antenna  system  both  have  infinitely  large  conduct- 
ivity, so  that  the  na  and  a are  both  infinitely  large  and  therefore, 
that  losses  are  negligible.  Also,  the  large  ground  plane  covered  with 
aluminum  foil  to  simulate  the  natural  ground  at  the  higher  frequencies 
is  assumed  to  influence  the  radiation  pattern  only  to  a minor  extent. 

Although  the  scale  models  are  built  mainly  for  radiation 
pattern  studies,  they  will  give  qualitative  information  on  the  impedance 
as  a function  of  frequency.  The  impedance  measurements  are,  there- 
fore, carried  out  in  the  laboratory  on  the  scale  model  between  200 
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MHz  and  1200  MHz  using  a HP  805A  slotted  line  and  a HP  415D 
VSWR  meter  for  frequencies  above  500  MHz,  and  a HP  803 A VHF 
impedance  bridge  with  a HP  417A  VHF  detector  for  the  frequencies 
below  500  MHz.  A difficulty  of  measuring  an  impedance  Z^  of  an 
unknown  load  when  ZT  is  much  greater  than  the  characteristic 
impedance  of  a transmission  line  Zq  was  pointed  out  in  Section  2. 2. 
As  one  way  to  eliminate  this  difficult} , the  line  length  d from  the 
position  of  the  load  Z L to  the  measuring  point  Z^  was  shortened. 
However,  on  the  scale  model  simulating  an  infinitely  large  natural 
ground,  the  measurement  equipment  could  not  be  mounted  right  at 
the  input  terminal  of  the  antenna,  thereby  -tracing  d to  zero. 

A 50- ohm  solid  copper  outer  conductor  coaxial  line  of 
approximately  20  centimeters  was  used  to  connect  the  antenna  to  the 
measuring  equipment.  Figure  2.  8 shows  a schematic  diagram  of  the 
experimental  arrangement.  Each  measurement  was  performed, 
first  with  a short  at  the  position  of  the  load  and  then  with  the  load 
connected.  This  technique  accounts  for  the  length  of  line  used  between 
the  measurement  equipment  and  the  load,  when  a Smith  chart  is  used 
to  obtain  the  actual  input  impedance  of  an  antenna.  The  measurement 
performed  in  this  manner  is  reliable  at  higher  frequencies  where 
Zl  = Zq  than  it  is  at  the  lower  frequencies  where  > > Z^. 
The  results  are  shown  in  Figs.  2. 15  and  2. 16. 

In  a further  study,  another  set  of  scale  models  were 
built  using  hemispherical  ground  planes,  instead  of  flat  ciicular 
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Simulating 

infinite 


Fig.  2.  8,  Block  diagram  showing  an  impedance  measure- 
ment set-up  for  a scale  model  antenna  system 
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ground  planes,  to  examine  any  similarities  in  the  radiation  patterns 
of  the  two  systems.  The  frequencies  selected  were  such  that  the 
diameter  of  the  hemispnerical  ground  planes  was  less  than  a wave- 
length at  1476  MHz.  This  is,  in  effect,  a 49. 2:1  scale  model  for 
simulating  operation  between  5 MHz  and  30  MHz.  Impedance  measure- 
ments were  also  made  on  this  model,  as  well  as  the  radiation  pattern 
studies  reported  in  Section  4. 4, 

In  studying  the  impedance  of  the  scale  models,  the  flat 
circular  disc  ground  planes  were  placed  at  heights  a^  and  a.^  above 
the  simulated  natural  ground.  The  displacement  is  the  radius 
of  the  semi-spherical  ground  plane.  The  displacement  is  the 
radius  of  the  flat  disc  equal  in  surface  area  to  that  of  the  spherical 
ground  plane. 

Figure  2.  9 shows,  in  detail,  the  experimental  arrange- 
ments of  scale  models.  In  the  first  case,  Fig.  2.  9(a),  the  displace- 
ment of  a finite  ground  plane  above  an  infinite  ground  is  the  same  as 
the  radius  of  a hemispherical  ground  plane  a.  In  the  second  case, 

Fig.  2.  9(b),  the  displacement,  a , of  a finite  ground  plane  above  an 
infinite  ground  is  chosen  to  be  equal  to  0.  707  times  the  radius  of  a 
circular  disc.  With  this  value,  the  surface  area  of  the  hemisphere, 
27ra2  , is  equal  to  the  surface  area  of  the  circular  disc.  That  is, 
in  the  first  case,  displacements  above  an  infinite  ground  plane  are 
equal,  and  in  the  second  case  the  surface  areas  of  the  small  ground 
equal.  The  object  of  this  study  was  to  determine  the  relative 
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impedance  characteristics  of  the  two  antenna  systems  and  to  determine 
the  dependence  of  the  system  upon  ground  plane  surface  area  and  upon 
its  location  with  respect  to  the  larger  ground  plane.  The  results  of 
the  laboratory  impedance  measurement  on  both  scale  models  are  shown 
in  Figs.  2. 15  and  2. 16. 

2. 4 Results  of  Measurement 

Antenna  input  impedance  measurements  taken  with  5,  2.  5, 

1. 25,  and  0. 625  meter  diameter  ground  plane  at  heights. 5,  2.  5,  1. 25, 
and  0. 625  meter  above  natural  ground  are  shown  as  a function  of 
frequency  in  Figs.  2. 11  to  2. 14.  The  data  are  arranged  to  show  the 
effect  of  ground  plane  sizes  with  their  location  fixed  at  constant  level 
and  the  effect  of  ground  plane  locations  with  the  size  fixed. 

The  capacitive  reactance  components  are,  in  general, 
less  dependent  upon  the  ground  plane  height  than  upon  its  diameter. 
However,  the  resistive  component  varies  with  the  height,  markedly 
and  shows  a peak  at  particular  frequencies.  The  input  impedance  of 
a monopole  antenna  that  is  located  over  an  infinite  ground  plane  is 
shown  in  Figs.  2. 10(a)  and  2. 10(b).  These  curves  are  derived  from 
Brown  and  Woodward’s  (Ref.  2)  experimental  data  where  the  length 
to  diameter  ratio  of  approximately  400  is  used.  Comparison  between 
results  of  infinite  ground  planes  and  finite  ground  planes  above  natural 
ground  shows  the  two  are  similiar  in  reactance  but  not  in  resistance. 

The  results  from  measurements  using  actual  size  of 
antenna  show  that,  for  a fixed  diameter  disc  ground  plane,  a resistive 


(a) 


Semi- spherical  surface  area  Ag  = 2ira2 

Circular  disc  surface  area  A~  = 2ir  a2  = A 

D s 


(b) 

Fig.2<9.  Scale  models  used  for  impedance  measurements 
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Frequency  (MHz) 

Fig.  2.  10(a). Input  resist?nce  as  a function  of  frequency 
for  a quarter  wavelength  monopole  at  30  MHz 
on  an  infinite  ground  plane 
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Fig.  2. 10  (b).  Input  reactance  as  a function  of  frequency 
for  a quarter  wavelength  monopole  at  30  MHz 
on  an  infinite  ground  plane. 
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Fig.  2.  16.  Input  impedance  versus  frequency 
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peak  occurs  at  higher  frequencies  as  the  location  of  this  ground  plane 
gets  higher  with  respect  to  an  infinite  ground.  Also,  the  peak  values 
of  the  resistive  components  are  decreasing  as  the  fixed  diameter  size 
is  increased.  A general  conclusion  which  can  be  drawn  from  these 
results  is  that  the  resonant  peak  occurs  at  a frequency  where 
k(-^-  + a)  is  constant.  Therefore,  as  a diameter  of  the  ground  plane 
and  a displacement  above  infinite  ground  plane  gets  larger,  a peaking 
occurs  at  lower  frequencies.  As  lor  the  reactive  components  of  the 
impedance,  the  reactance  versus  frequency  curve  is  similar  to  that 
of  a monopole  above  an  infinite  ground  plane  when  the  diameter  of  a 
disc  is  5 meters.  The  displacement  above  an  infinite  ground  does  not 
seem  to  affect  the  reactive  component  as  much  as  the  diameter  of  a 
finite  ground  plane.  However,  even  for  D = 5 meters,  reactance 
values  are  lower  than  that  of  a monopole  antenna  with  the  same  length 
to  diameter  ratio  on  an  infinite  ground. 

A limited  number  of  impedance  measurements  has  been 
conducted  using  a scale  model  which  is  mainly  constructed  to  take  a 
far-zone  radiation  pattern  of  the  antenna  and  ground  plane  system. 

The  input  terminals  as  well  as  infinite  natural  ground  simulations 
have  not  been  able  to  scale  ideally.  These  discrepancies  usually 
affect  impedance  measurements  more  than  radiation  patterns  obtained 
at  far-zone  area.  Under  these  circumstances,  only  a qualitative 
comparison  of  the  impedance  has  been  possible. 

First,  the  actual  model  and  the  scale  model  antenna  and 
ground  plane  system  at  a ground  plane  diameter  equal  to  a half  wave- 
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length  at  30  MHz  compared  favorably  in  both  real  part  and  imaginary 
part. 

Secondly,  the  scale  model  antenna  with  a disc  ground  plane 
and  a hemispherical  ground  plane  comparison  shows  that  neither  the 
diameter  of  the  two  types  of  ground  planes  or  the  area  of  the  ground 
planes  has  a basis  of  similarity  in  impedances.  However,  it  is 
deduced  qualitatively  from  the  limited  number  of  experimental  data 
that  the  distance  between  the  base  of  the  antenna  to  the  natural  ground 
and  the  conductivity  along  this  path  has  a greater  bearing  upon  the 
resonance  phenomena  observed  in  the  impedance  measurements.  In 
other  words,  k(a  + y-)  = const  for  a disc  ground  plane  and  ka  = 
const  for  a semi- spherical  ground  plane  will  determine  the  resonance 
conditions.  These  constants  seem  duferent  in  general  for  different 
ratio  of  a to  D. 

In  addition  to  these  observations,  the  determination  of 
whether  resonance  peak  observed  in  input  resistance  measurement  is 
largely  due  to  increase  in  radiation  resistance  or  not  will  be  shown 
in  the  rest  of  this  chapter  and  in  Chapter  5. 

In  studying  impedance  measurement  data,  it  is  also 
observed  that  both  real  and  imaginary  part  of  the  impedance  behave, 
in  most  cases,  such  as  in  Fig.  2. 17(c).  Frequencies  where  peaking 
effect  occurs  are,  of  course,  a function  of  a and  D.  Figure  2. 17(a) 
shows  an  input  impedance  as  a function  of  frequency  for  a monopole 
on  an  infinite  conducting  ground  plane  and  its  equivalent  circuit. 
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Figure  2. 17{b)  shows  a parallel  resonance  circuit  whose  impedance 
characteristics  add  up  with  Fig.  217(a)  to  give  an  impedance  char- 
acteristic which  is  observed  in  the  measurements. 

Therefore,  it  is  possible  to  synthesize  an  equivalent 
circuit  of  a type  shown  in  Fig.  2. 17(c)  to  further  study  the  effect  of 
a finite  ground  plane  upon  an  equivalent  circuit.  In  this  way,  a cor- 
relation may  be  obtained  between  a,  D and  the  circuit  parameters. 

2.  5 Copper  Losses  Due  to  the  Antenna  and  the  Ground  Plane 

The  real  part  of  the  antenna  input  impedance  contains  a 

part 

Rr  = “T  & J Re  (E  x H*)  • dS  (2.9) 

that  is  directly  proportional  to  the  radiated  power  of  an  antenna,  for  a 
constant  input  current.  Consequently,  it  is  important  to  separate 
the  resistive  component  of  the  input  impedance  measurement  data 
into  the  resistive  loss,  which  is  due  to  several  causes,  as  explained 
in  the  following  section,  and  the  radiation  resistance.  This  permits 
determination  of  whether  the  unusual  variation  of  the  input  resistance 
as  a function  of  frequency,  which  has  been  found  in  the  input  impedance 
measurements,  is  due  to  an  increase  in  loss  or  to  an  increase  in  the 
radiation  resistances  at  particular  frequencies. 

The  total  antenna  resistance  is  the  sum  of  the  several 


separate  components 
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(1)  Radiation  resistance  Rr 

(2)  Ground  terminal  resistance  R 

g 

(3)  Resistance  of  tuning  units  R^ 

(4)  Resistance  of  equivalent  insulation  loss  R. 

(5)  Resistance  of  equivalent  conductor  losses  R^, 

(6)  Transmission  line  loss  R . 

Among  the  six  separate  contributions,  R^  and  are 

not  present  in  this  case  because  there  are  no  tuning  units  attached  when 

the  measurements  are  taken  and,  as  a low  potential  receiving  antenna, 

no  appreciable  insulation  losses  are  involved. 

The  loss  due  to  the  transmission  line  (R  ) was  evaded 

m 

by  measuring  the  impedance  at  the  base  directly  below  the  ground 
plane,  for  the  actual  size  and  by  taking  a short  circuit  measurement 
with  the  shcrt  placed  at  the  input  terminal  of  the  scale  model.  The 
loss  in  the  line  was  subtracted  from  the  measured  data. 

In  the  out-door  measurement  with  the  natural  ground 

below  the  finite  aluminum  ground  plane  R cannot  be  exactly  calcu- 

g 

lated  without  precise  information  of  the  conductivity  and  other  para- 
meters of  the  dirt  ground.  The  impedance  measurement  on  the  scale 
model  when  the  natural  ground  is  simulated  by  the  aluminum  foil 

enables  both  R and  R to  be  computed, 

gw 

The  following  analysis  is  to  permit  the  calculation  of  R^ 

for  both  the  actual  model  and  the  scale  model  and  R of  the  scale 

g 


model. 


41 


2.  5. 1 The  Internal  Impedance  of  the  Plane  Conductor. 

The  current  density  resulting  from  the  movement  of  charges  in  a 
conductor  is  given  by  Ohm’s  law: 

J = oE  (2.8) 

The  constant  a is  the  conductivity  of  the  conductor,  and  the  Maxwell's 
equation  is 

VxH  = aE  + f-  (2.9) 

dt 

For  a harmonically  oscillating  field  with  time  dependence 

V x H = (a  + juje)E  (2. 10) 

In  the  absence  of  free  charges  p,  V • D = 0.  Also,  for  most  con- 
ductors, the  displacement  current  3D  / 3 1 is  negligibly  small  com- 
pared to  the  conduction  current. 

Then, 

VxVxE  = V(V  • E ) - V2  E = Vx  (-1*1)  = JL  n V X H 

(2.11) 

•*.  VaE  = = ja>  paE 

Similarly,  using  J = oE 

V2  J = jcupoj 


(2.  12) 
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^ i Conducting  Surface 


/ 


z 


Fig.  2. 18.  Plane  solid  conductor 

For  a plane  conductor  of  infinite  depth,  with  no  field  vari- 
ations along  the  length  or  width  as  shown  in  Fig.  2. 18. 


d2  i 

— = ju>poi 
dx2  z 


(2. 13) 


where 


T2  = jWpO 


Since 


ILLil 

yfT 


and 


t = (1  + j)  yjTifi'j.0  (2.14) 

where  by  definition  6 - — - — , is  called  depth  of  penetration 

^ Tll^O 

of  the  field,  or  the  skin  depth. 

The  solution  of  differential  equation  2.  13  is  then  given  as 
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l = 


*0 


-d  f j). 
6 ‘ 


(2. 15) 


The  total  current  flowing  in  the  plane  conductor  is  found  by  integrating 
the  current  density  i from  the  surface  to  the  infinite  depth.  For  a 
unit  width, 

(l  + j)  . 

°o  - x in6 

i dx  = f in  e dx  = -= — — (2. 16) 

z ^0  1 + ] 

The  electric  field  at  the  surface  is  E = . Therefore,  the 

zo  o ’ 

internal  impedance  per  unit  length  and  unit  width  is 


1 + j 
ob 


(2.  17) 


If  Z is  defined  as  Z = R + jtuL.,  then 
s s s i’ 


(2.  18) 


Rg  is  the  resistance  of  the  plane  conductor  for  a unit 

length  and  unit  width.  For  a finite  area  of  conductor,  the  resistance  is 

obtained  by  multiplying  R by  length,  and  dividing  by  the  width. 

s 

For  a circular  ground  plane  with  radial  current  distribution, 

the  total  surface  resistance  is  obtained  by  multiplying  R by  the 

s 

radius  and  dividing  by  the  mean  circumference  of  the  plate. 

For  an  aluminum  ground  plane  which  has  material  constants 


of: 
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rt 

a = 3.72  x 10  mhos/ meter 

n 

p = 4ir  x 10  henries/ meter 

* 0.0826 

6 = zz  meters 

VT 

f = frequency  in  Herta 
the  surface  resistivity  is  computed  to  be 

R = 3.26  x 10"7  VT  (2. 19) 

s 


which  becomes 
at  f = 30  MHz 

= 5 MHz 


Rg  = 3. 26  x 10~7  a x 10 7 = 1.  7g_xU0lj-ohms/ 

square 


= 3.  26x  10" 7 ^5  x 106  = 7.3  xlO'4 


When  the  aluminum  ground  plane  has  a diameter  equal  to  a 
quarter- wavelength  at  30  MHz,  the  radius  r is  2.  5 meters  and  the 
mean  circumference  is  Ttr.  The  total  surface  resistance  is,  therefore, 


R 

U) 


(2.20) 


The  numerical  values  computed  at  5 MHz  and  30  MHz 


become 
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R , = — ohms  = 5.  7 x 10  ^ ohms  at  30  MHz 

Cl)  O.  14 

-4 

7 3 x 10  -4 

= - g — ohms  = 2.32  x 10  ohms  at  5 MHz  . 

The  contribution  of  the  ground  plane  surface  resistance  toward  the 
input  resistance  is,  therefore,  negligible. 


2.  5. 2 Internal  Impedance  of  a Conductor  with  a Circular 


Cross  Section 


2.  5.2. 1 Current  in  a Wire  of  a Circular  Cross  Section. 


Let  the  current  flowing  on  the  antenna  of  a circular  cross  section, 


1/4  inch  in  diameter  be  assumed  to  flow  mainly  on  the  axial  direction; 

F = i 'z . Also,  no  axial  or  circumferential  variation  is  assumed, 
z ' 


Fig.  2. 19.  Current  in  a cylindrical  wire 
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Then  V2  i = jcuptai  becomes  in  the  cylindrical  coordinate  system, 


d2i  i di 

-~rV  + ~ — r2-  = iu>/ioi 

dr2  r dr  J z 


(2.21) 


If  we  let  T2  = -jto/icr 


d2i 

ZL 

d 2 
z 


1 di2 


r dr 


T2  i 


= 0 


(2.22) 


For  a solid  wire,  the  solution  must  be  finite  at  r = 0 . 
Therefore,  it  takes  the  form  of 


*z  = AJo(Tr) 


(2.23) 


Let 


lz  = lo  at  r = r0 


(2.  24) 
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then 


A = 


‘0 

w 


(2.  25) 


1 = 1 


J0(Tr> 

0 w 


(2.  26) 


Since 


rp  2 

T = -JO Olio 


(2.27) 


Since 


Ber  (v)  f j Bei  (v) 


(2.28) 


Where 


Ber  (v)  h Real  part  of  JQ  (vt) 
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Bei  (v)  s Imaginary  part  of  J-V  v \ 

V ) 


The  current  density  in  the  axial  direction  can  be  written 


as  now 


Ber^T~)  * > Bei  (^T~) 

Be'\ir-) + ) 


(2.  29) 


If  the  ratio  of  rQ/6  is  large,  the  I i^ / iQ  plot  will  agree 


closely  with  the  plane  conductor  derivation  of 


X-  e'(r0'r>/6 


(2.30) 


where  r^  - r replaced  x for  the  case  of  a plane  conductor. 
Also, 


$ H • df  = I and  2?rrnH^  = I 

y 00  r=rQ 


(2.31) 


From  the  Maxwell's  equations, 


V x E - -jcu/iH 


(2.32) 


and  for  the  round  wire 


i dE 

l - — 1 — — JL 
0 jw  M dr 


(2.33) 
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Since 


‘o  J0(Tr> 

' TO 


Htf  = - 


«0  J0 (Tr> 

T JgiTr^5 


Consequently, 


I 


2ttt 


_0  t J(J  <Tj*> 

T 0 ~j^[ T?A) 


(2.34) 


(2.35) 


(2.36) 


2.  5.  2.  2.  The  Internal  Impedance  of  a Round  Wire. 
The  internal  impedance  Z i is  defined  as 


„ . Ezlr=ro  TJo(Tro> 

i - I = 2”VJ0(Tr0> 

Using  the  formula  that 


Ber  (v)  + j Bei  (v) 


(2.37) 


and 


Ber’(v)  + j Bei'(v) 


--  [Ber  (v)  + j Bei  (v)J 


(2.  38) 
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The  internal  impedance  becomes 


Z.  = R + itoL.  = 

1 1 y[2u 


r Ber  q + j Bei  q 1 (t> 

LBer’ q+ j Bei’qJ  ( * y) 


where 


■ TT 


(2.  40) 


(2.  41) 


Ber  q Dei*  q - Bei  q Ber1  q ohms/  meter 
(Ber'  q)2  + (Bei  q)2 


Ber  q Ber’  q + 
(Ber'  q)2  + 


(2.  42) 


Bei  q Bei1 
(Bei  q)2 


ohms/ meter 


Using  the  same  analysis  except  that  the  wall  thickness  is 
small  enough  with  respect  to  the  radius  of  the  tube  to  be  able  to  consider 
the  tublar  conductor  as  a flat  conductor  of  finite  thickness,  the  internal 
impedance  can  be.  found  to  be 


Z = (1  + j)  R, 


Cosh  rd 
Sinh  rd 


(2.43) 


j 
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where 


T = Sl+JL 

a 


d = thickness  of  the  tubular  wall 

r = ViSL 

s a 


Ther 


ohm/ unit  length 
(2.44) 


The  impedance  per  unit  length  is,  then. 


R 


“Li  --2^ 


0 


Sinh  \2d/Y 


Sin  \26^  5 1 


Cosh  rd/6 


cos  Kai 


R 


R = 


2ur 


0 


Sinh  (2d^  s)  + Sin  (2d/  fi) 
Cosh^d  sj  - Cos^d  sJJ 


ohm/ unit  length 


(2.  45) 


ohm/ unit  length 


The  tubular  antenna  used  for  impedance  measurement  has 
the  following  dimensions: 

_ 3 

radius  = 3.17  x 10  meters 

-4 

wall  thickness  - 6.31x10  meters 
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length  of  antenna  = 2.  5 meters 
The  skin  depth  6 computed  for  the  copper  at  both  ends  of 
the  frequency  band  are 

6 = 1. 1 x 10  ^ meters  at  30  MHz 


= 3 x 10  ^ meters  at  5 MHz 
Therefore,  the  ratio  2d/ 6 for  equation  2.45  becomes 

2d. '6  = 11.4  at  30  MHz 
= 42.4  at  5 MHz 

For  large  values  of  x , both  sinh  x and  cosh  x approach 
(1/  2)e*  and  sinh  x ' , sin  x and  cosh  x ' „ cos  x . Equation 
2.45  then  simplifies  to 

R . / 

R ~ coL.  when  ''  6 >1  (2.46) 

w 27ir„  l 


For  the  copper  tubular  antenna  2.  5 meters  long  the  real 
part  of  the  internal  impedance  becomes 


R = 0.  189  ohms  at  30  MHz 

U) 


: 0.  063  ohms  at  5 MHz 

The  contribution  of  the  copper  loss,  as  shown  above,  toward 
the  antenna  input  impedance  is  also  nog’ -gillie. 


53 


In  this  section,  a numerical  calculation  of  ground  terminal 

resistance  R and  conductor  loss  resistance  R was  made.  In  the 
g w 

act’ial  antenna  and  ground  plane  system,  the  conductivity  of  the  infinite 

ground  beiow  the  finite  ground  must  be  known.  However,  calculations 

pex  formed  here  did  net  lake  the  natural  ground  into  account.  The  only 

ground  considered  is  an  aluminum  ground  plane  used  as  a finite  ground 

below  the  antenna.  Contributions  from  the  tubular  antenna  of  2.  5 

meters  in  length,  used  for  a monopole,  to  R^  were  also  considered. 

In  the  frequency  range  of  5 to  30  MHz,  where  the  experiment  was 

carried  out,  the  ground  loss  from  the  largest  ground  plane  (D  = 5 

-4 

meters)  used  was  computed  to  be  less  than  5.  7 x 10  ohms.  The  loss 

from  the  tubular  antenna  conductor  itself  was  less  than  0.  2 ohms. 

In  the  scale  model,  the  ground  plane  used  to  simulate  an 

infinite  natural  ground  was  a 4-foot  diameter  aluminum  circular 

ground  plane.  The  conductivity  data  and  skin- depth  data  available  for 

aluminum,  both  large  and  small  ground  plane  loss  were  calculated. 

Also,  a conduction  loss  from  an  AWG  No.  28  wire  of 

6.  25  cm  in  length  used  to  simulate  a monopole  antenna  was  calculated. 

The  overall  loss  due  to  R and  R for  the  scale  model  was  well 

g w 

below  1 ohm  and  there  was  no  obvious  peaking  effect  due  to  R and 

R^  as  frequency  was  varied  from  200  MHz  to  1200  MHz. 

Therefore,  it  seems  that  the  resonant  peak  observed  in  the 

real  part  of  the  antenna  input  impedance  in  the  experiment  originate 

from  sources  other  than  R and  R . 
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CHAPTER  III 


CURRENT  MEASUREMENTS 


3.  1 Introduction 

In  this  chapter,  an  extensive  study  of  the  current  distribu- 
tion on  a monopole  antenna  one  quarter  wavelength  long  at  30  MHz  is 
reported  as  functions  of  a ground  plane  size  and  its  location.  King  and 
Harrison  (Ref.  9)  computed  analytically  the  current  distribution  of 
a symmetrically  center- driven  antenna.  For  this  dipole,  the  current 
distribution  was  found  to  be  extremely  close  to  the  sinusoidal  distri- 
bution that  is  normally  assumed,  so  long  as  the  antenna  length  to  di- 
ameter ratio  is  very  large.  However,  the  effect  of  a small  ground 
plane  whose  diameter  is  less  than  a half  wavelength  upon  the  antenna 
current  when  the  system  is  placed  above  an  infinite  ground  has  not 
been  reported. 

The  particular  emphasis  here  is  to  show  the  differences 
between  the  current  distribution  on  an  antenna  over  an  infinite  ground 
plane  and  the  same  size  of  antenna  over  a finite  ground  displaced  less 
than  a wavelength  above  an  infinite  ground  plane. 

3.  2 Theory  of  Current  Probe 

In  order  to  detect  the  current  densities  on  the  antenna,  a 
rectangular  loop  antenna  whoso  dimensions  are  much  smaller  than  a 
wavelength  was  used  as  a probe.  The  theory  of  operation  of  this  loop 
as  a probe  has  been  studied  by  King  and  Whiteside  (Ref.  10  ).  The 


54 


CHAPTER  HI 


CURRENT  MEASUREMENTS 


3. 1 Introduction 

In  this  chapter,  an  extensive  study  of  the  current  distribu- 
tion on  a monopole  antenna  one  quarter  wavelength  long  at  30  MHz  is 
reported  as  functions  of  a ground  plane  size  and  its  location.  King  and 
Harrison  (Ref.  9)  computed  analytically  the  current  distribution  of 
a symmetrically  center-driven  antenna.  For  this  dipole,  the  current 
distribution  was  found  to  be  extremely  close  to  the  sinusoidal  distri- 
bution that  is  normally  assumed,  so  long  as  the  antenna  length  to  di- 
ameter ratio  is  very  large.  However,  the  effect  of  a small  ground 
plane  whose  diameter  is  less  than  a half  wavelength  upon  the  antenna 
current  when  the  system  is  placed  above  an  infinite  ground  has  not 
been  reported. 

The  particular  emphasis  here  is  to  show  the  differences 
between  the  current  distribution  on  an  antenna  over  an  infinite  ground 
plane  and  the  same  size  of  antenna  over  a finite  ground  displaced  less 
than  a wavelength  above  an  infinite  ground  plane. 

3.  2 Theory  of  Current  Probe 

In  order  to  detect  the  current  densities  on  the  antenna,  a 
rectangular  loop  antenna  whoso  dimensions  are  much  smaller  than  a 
wavelength  was  used  as  a probe.  The  theory  of  operation  of  this  loop 
as  a probe  has  been  studied  by  King  and  Whiteside  (Ref.  10  ).  The 


54 


56 


components  are  also  preb.mt,  even  though  small,  which  are  proportional 
to  the  average  electric  field  in  the  plane  of  the  loop. 

When  electro-  magnetic  fields  are  incident  upon  a loop,  there 
are  re- radiated  fields  in  addition  to  the  incident  field.  The  total  field, 
therefore,  is 

"B  = B1  + Br  (3.1) 

where  ¥ is  the  total  field,  "B1  the  incident  field,  and  B*  the  re-ra- 
diated  field.  From  Maxwell’s  curl  equation, 

V x E - - ^ (3. 2) 

iU)t 

Assuming  that  the  incoming  signal  has  eJ  time  dependence,  the 
curl  equation  becomes 

VxE  = -ju?B  - -jkcB  (3.3) 


where 


k 


277 

\ 


c 


Integrating  both  sides  over  loop,  this  becomes 

ff  VxE  - dS  = - jk  Jf  c B • dS  (3.  4) 

(3.  5) 


(6  E • df  = -jkc  Jf  B ’ dS 
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If  the  loop  area  is  very  small,  then  the  magnetic  field  in 
the  loop  is  uniform  with  magnitude  B1^  . The  integral  is 

e = -jkc  fj  3 • d§>  = -jkc  Bq  A (3.6) 

Where  e is  the  induced  emf;  and  where  A is  the  area  of  a loop.  De- 
fining an  effective  height  of  a loop  antenna  as 

i 

1^  = -jkA  (3.7) 

the  low  frequency  input  admittance  of  a loop  with  a constant  current 
(Ref.  12)  is 

Y0  = [^v1  +i$r  ft  -ir)kR  d7'  ■ d7]  1 (3-8) 

An  unloaded  magnetic  sensitivity  Kg  is  defined  as 

KB  = YoVA  (Ref-  10 } (3.9) 

Now 

e = y Yq  = -jkAcB0  (3.10) 

I0  - (-JkA)  Yq  cB0  - Y0  CB0  . AKb(cB0)  (3.  11) 

IQ  , therefore,  is  a measure  of  the  normal  component  of  the  magnetic 
field  at  the  center  of  the  loop,  independently  of  the  electric  field.  The 
loop  sensitivity  for  the  electric  field,  Kg  , is  a function  of  the  size 
of  the  loop. 
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Let  the  current  in  the  loop  due  to  the  incident  fields  consist 
of  the  zero  and  the  first  phase- sequence  currents,  i.  e. , 


I(s)  = I°(s)  + I^(s) 

I(s  +£/2)  = I°(s)-  I(1)(s) 


(3. 12) 


y 

i 

0(S) 

1/2 

S = ( 

Fig.,  3.  1 Rectangular  current  loop 

The  first  phase  current  1^  (s)  depends  directly  upon  the  electric  field 
(Ref,  10).  Let  the  loop  lie  on  the  x-y  plane  and  the  center  of  the  loop 
at  the  origin,  then  (s)  can  be  further  broken  up  into 

I*1)  = I(1^  + I*1*  (3.  13) 

x y 

corresponding  to  E ^ and  E ^ . Therefore,  the  current  at  s = 0 

X0  y0 
is 

I(1\o)  = h Y.  E 1 (4) 

y * i y0 

n 

For  a short  dipole,  h^  = /2  where  £ is  the  length  of  the  antenna. 


(3. 15) 


r*'(0)  = X K-,  E ' ' 
y E yQ 


where 


KE  - YIhel/X 


(3. 16) 


For  a magnetic  loop  with  a diameter  D , the  effective  height  hfa  becomes: 


1^  = -jku  D2/  4 


(3. 17) 


and  the  magnetic  sensitivity  defined  in  Eq.  3.  9 becomes 


kb  ■ 


7rD 


At iJa  -3.52+  13.0—) 

°\  X2/ 


(3. 18) 


where 


fl  = 2 In 


2D 


a = radius  of  wire  loop 

The  current  induced  on  the  portion  of  the  loop  parallel  to 
the  incident  field  is  from  the  electric  sensitivity. 


1(0)  = E 


i D 

y0  mx 


COS  0 


(3.19) 


i 
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where  tj  is  the  characteristic  impedance  of  the  medium,  and  a^  is 
an  expansion  coefficient  given  by  Storer  (Ref.  13).  For  D 0.  IX  , 
the  electric  sensitivity  reduces  to 

K£  = (j  2jt2  D2/A  2)/  770  ( n-  3.  52)  ( 1 - 9.8D2/X  *)  (3.20) 


Fig.  3.  2.  Circular  loop  probe 


Also,  it  can  be  shown  that  for  a loop  with  a diameter  D < < 0.  03  X 


Ke  oc 


Kb  oc 


(3.21) 


From  the  dependences  of  electric  and  magnetic  sensitivities 
upon  the  loop  size,  and  since  the  loop  current  is  directly  proportional 
to  the  sensitivities  Kg  and  Kg  , the  loop  diameter  should  be  as  small 
as  possible  to  minimize  the  effect  of  electric  field  upon  the  total  current 
measured.  The  lower  limit,  however,  is  dictated  by  the  fabrication 
technique  and  the  receiver  sensitivity. 
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3.  3 Experimental  Procedures 

The  current  measurement  on  the  antenna  was  obtained  by 
detecting  the  magnetic  field  locally  along  the  antenna  with  a circular 
probe  that  is  sensitive  to  magnetic  fields.  Initially,  the  experiment 
was  carried  out  on  a model  at  a frequency  range  between  5 MHz  through 
30  MHz.  The  antenna  used  for  the  impedance  measurement  was  2.  5 
meters  long,  which  represents  a quarter  wavelength  at  30  MHz.  A 
one- quarter- inch  copper  tubing  was  used  an  an  antenna,  and  the  current 
along  the  antenna  was  detected  with  a traveling  probe  moved  along  an 
axial  slit  cut  in  the  antenna.  This  technique  was  first  developed  by 
King  et  abfor  measurement  of  antenna  current  distributions. 

Since  the  object  of  this  study  was  to  examine  the  effects 
of  the  size  and  location  of  a finite- size  ground  plane,  a systematic 
variation  of  ground  plane  diameters  and  their  locations  above  the 
natural  ground,  up  to  a half  wavelength  at  the  highest  frequency,  was 
planned.  In  the  frequency  band  chosen,  this  size  represented  a ground 
plane  diameter  and  its  location  up  to  5 meters.  Due  to  the  bulky  physi- 
cal size  of  the  full  scale  system,  a smaller  scale  model  wac  adopted. 

The  scale  models  used  were  the  same  ones  used  in  the 
radiation  pattern  studies.  For  current  distribution  measurements, 
the  entire  system  was  put  into  an  anechoic  chamber  that  is  designed 
to  operate  above  100  MHz.  The  detecting  probe  for  the  local  magnetic 
field  along  the  antenna  has  a loop  diameter  much  smaller  than  a wave- 
length, even  at  the  nigher  frequency  of  200  MHz  through  1200  MHz. 
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The  probe  was  suspended  from  the  top  of  the  anechoic  chamber  to  the 
current  element  and  its  movement  was  synchronized  using  a Scientific- 
Atlanta  recorder.  The  detected  signal  was  then  fed  into  a recorder  and 
the  result  recorded. 

The  probe,  the  anechoic  chamber,  and  the  equipment  used 
for  this  measurement  are  shown  in  Fig.  3.  3 through  3.  5.  The  block 
diagram  showing  the  current  measurement  set-up  is  shown  in  Fig.  3.  6. 

Detection  of  a current  along  a current  element  with  an 
external  probe,  as  shown  in  our  experiment,  was  accepted  fully  by 
many  people.  Since  the  near  field  may  be  disturbed  by  the  presence 
of  the  probe,  and  its  lead,  a probe  traveling  inside  an  antenna  with 
only  that  part  of  the  probe  which  is  used  for  detecting  the  local  mag- 
netic field  exposed  through  a slit  was  the  optimum  choice.  However, 
if  the  probe  is  very  small,  compared  with  a wavelength,  and  if  the 
lead  from  the  probe  to  the  recorder  is  well  shielded  from  any  stray 
pick-up,  the  effect  of  the  external  probe  upon  the  antenna  current  was 
minimal.  The  lead  from  the  probe  to  the  recorder  was  a solid  copper 
outer- conductor  coaxial  line  (Micro  Coax  cable)  with  an  outer  diameter 
of  0.  5 mm.  The  probe  lead  was  oriented  perpendicular  to  the  antenna 
current  element,  to  minimize  the  disturbance  of  the  near- zone  electric 
field  parallel  to  the  current  element. 

To  test  the  accuracy  of  this  type  of  current  measurement, 
a quarter- wavelength  monopole  was  mounted  over  a ground  plane,  that 
simulated  natural,  infinite  ground.  The  ground  plane  was  a 5 wavelength 
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Scientific- Atlanta 
recorder 


Fig.  3. 6.  Block  diagram  for  current  measurement 
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in  diameter  aluminum  ground  plane.  The  test  frequency  was  1200  MHz, 
where  6. 25  cm  is  a quarter- wavelength.  The  measured  current  dis- 
tribution at  this  frequency  was  plotted  against  an  exact  sinusoidal  dis- 
tribution in  Fig.  3.  7.  This  measurement  agrees  very  closely  with  the 
King  measurement  of  a current  on  a cylindrical  antenna  with  a probe 
traveling  inside  the  cylindrical  antenna  with  a slit  cut  along  the  axial 
direction. 

3.  4 Measurement  Results 

The  variables  for  the  current  measurements  were  the  same 
as  those  for  the  radiation  pattern  measurement.  In  Fig.  3.  8,  a plot 
of  the  theoretical  current  distribution  of  a monopole  over  an  infinitely 
large  conducting  ground  is  made  for  6 different  frequencies.  The  an- 
tenna was  a quarter -wavelength  at  the  upper  end  cf  the  frequency, 

namely  30  MHz,  and  as  frequency  was  reduced  to  5 MHz  in  4 steps, 
the  resulting  antenna  current  distribution  was  plotted  in  order  to 
compare  these  currents  with  those  on  the  same  length  of  antenna 
over  a finite  ground  plane  located  various  distances  away  from  the 
large  ground.  In  Figs.  3.9  through  3.14  at  each  frequency  the 
diameter  of  a small  ground  plane  was  fixed  at  a given  fraction  of 
wavelength  and  locations  of  this  ground  plane  were  varied  in  four 
steps.  Figures  3.  15  through  3.  20  show  the  effect  of  ground 
plane  diameter  when  located  at  a fixed  distance  away  from  the  large 


g 


Current  Probe  Location  From  the  Antenna  Base 
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Fig.  3.7.  Comparison  of  a sine  curve  with  a current  distribution  of 
a monopole  over  a large  ground  plane 


Current  Probe  Location  From  the  Antenna  Base 
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The  same  measurement  was  performed  with  a hemispherical 
ground  plane,  in  order  to  study  the  validity  of  a sinusoidal  current  as- 
sumed on  the  antenna  for  a theoretical  study  which  was  carried  out  in 
Chapter  5.  These  measurement  results  are  plotted  in  Figs.  3. 21 
through  3. 22. 

Only  amplitudes  of  the  current  on  the  antenna  have  been 
measured,  and  the  maximum  amplitude  for  each  set  of  measurements 
was  normalized  to  1. 

The  theory  developed  in  Chapter  5 for  an  antenna  on  a hemi- 
spherical ground  plane  began  with  an  assumption  that  the  current  dis- 
tribution on  the  antenna  is  in  the  sinusoidal  form.  It  was  further  as- 
sumed that  the  antenna  current  was  independent  of  the  induced  current 
on  the  ground  plane  surface.  Accuracy  of  the  theoretical  results  largely 
depend  upon  the  accuracy  of  this  assumption.  From  the  measurement 
results  shown  in  Fig.  3. 21  through  3.  22,  it  can  be  seen  that  at  frequen- 
cies where  the  antenna  is  near  a quarter- wavelength  long,  the  assump- 
tion of  a sinusoidal  current  distribution  was  fairly  accurate.  However, 
at  frequencies  below  20  MHz,  the  currents  on  the  antenna  were  modi- 
fied appreciably  due  to  interaction  between  the  antenna  current  and  the 
ground  plane  current.  At  5 MHz,  in  particular,  the  current  distri- 
bution on  the  antenna  was  almost  constant  in  magnitude,  similar  to  the 
current  distribution  normally  assumed  for  a Hertzian  dipole.  Origi- 
nal assumption  of  a sinusoidal  current  distr  ibution  should,  on  the  other 
hand,  provide  a form  close  to  a triangular  distribution  at  5 MHz 


71 


for  an  antenna  resonating  at  30  MHz.  This  type  of  distribution  is  simi- 
lar to  the  current  distribution  of  so-called  Abraham  dipole.  Assuming 
that  the  physical  length  of  these  two  antennas  were  the  same,  their 
effective  heights  differ  by  a factor  of  2.  Therefore,  a constant  current 
distribution  gives  a radiation  resistance  as  large  as  4 times  that  of  a 
triangular  current  distribution  provided  that  the  amplitude  of  the  tri- 
angular currents  are  the  same. 

A monopole  with  a finite  disc  antenna  shows  a better  agree- 
ment between  the  assumed  sinusoidal  current  distribution  and  the  ex- 
perimentally measured  distribution  at  both  the  high  and  low  ends  of 
frequency  band.  However,  at  the  midband,  the  current  distributions 
show  lower  amplitudes  compared  with  a sinusoidal  distribution.  Also 
noticed  was  that  for  a given  size  of  disc  diameters,  the  changes  in  the 
location  of  the  ground  plane  height  (a)  affects  the  current  distribution 
more  than  the  changing  ground  plane  diameters  at  a fixed  position.  It 
was  also  noticed  that  when  the  ground  plane  diameter  was  less  than 
1.  25  meters,  the  current  distribution  was  such  that  the  amplitudes  were 
bigger  than  those  of  sinusoidal  distribution. 
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Fig.  3.  14.  Current  distribution  on  a monopole  antenna  with  a 
finite  ground  plane  at  various  locations  with  respect 
to  an  infinite  ground  at  7.  5 MHz 
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Fig.  3. 15.  Current  distribution  on  a monopole  antenna,  with  a 
ground  plane  of  various  diameters  at  a given 
location  with  respect  to  an  infinite  ground  at  30  MHz 
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Fig.  3. 17.  Current  distribution  on  a monopole  antenna  with  a 
ground  plane  of  various  diameters  at  a given  loca- 
tion with  respect  to  an  infinite  ground  at  20  MHz 
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Fig.  3.  18.  Current  distribution  on  a monopole  antenna  with  a 
ground  plane  of  various  diameters  at  a given  loca- 
tion with  respect  to  an  infinite  ground  at  15  MHz 


Fig.  3.  20.  Current  distribution  on  a monopole  antenna  with  a 
ground  plane  of  various  diameters  at  a given  loca- 
tion with  respect  to  an  infinite  ground  at  7.  5 MHz 
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CHAPTER  IV 

RADIATION  PATTERN  MEASUREMENTS 


4. 1 Introduction 

The  important  parameters  of  any  antenna  are  its  radiation 
patterns  and  its  input  impedance.  In  this  chapter,  the  experimental 
results  of  the  radiation  pattern  measurements  are  given.  The  arrange- 
ment of  the  experimental  results  is  intended  to  show  the  effect  of  the 
ground  plane  size  and  its  location  with  respect  to  the  infinitely  large 
natural  ground  upon  the  radiation  properties  of  the  antenna  system. 
Certain  approximations  and  assumptions  had  to  be  made  because  of  the 
practical  limitations  that  confront  this  type  of  experiment,  and  they  are 
discussed. 


4.  2 Measurement  Problems 

The  input  impedance  measurements  discussed  in  Chapter  II 
were  performed  with  a monopole  a quarter  wavelength  long  at  30  MHz 
and  a ground  plane  whose  diameter  varied  from  0. 625  meter  up  to  5 
meters.  The  experiment  was  conducted  in  the  frequency  band  of  5 MHz 
through  30  MHz.  In  performing  the  radiation  pattern  measurements, 
several  assumptions  and  approximations  were  made.  The  antenna  test 
range  was  originally  designed  for  high  frequency  operation,  mainly  be- 
cause the  required  physical  size  and  distance  between  the  transmitting 
and  receiving  antenna  were  more  practically  realizable  in  the  higher 
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frequency  ranges.  The  linear  dimensions  of  the  original  antenna  and 
ground  plane  were  reduced  by  a factor  of  40  so  that  the  operating  test 
frequency  was  between  200  MHz  to  1200  MHz.  The  electromagnetic 
scaling  and  conditions  under  which  this  is  valid  have  been  discussed  in 
detail  in  Chapter  2. 

In  the  original  antenna  and  ground  plane  measurement,  the 
infinite  ground  below  a finite  ground  plane  was  the  natural  ground  with 
a certain  electrical  conductivity  (a  ).  One  of  the  conditions  for  scaling 
the  electromagnetic  system  is  increasing  the  conductivity  of  the  scale 
model  by  a scale  factor.  In  this  case,  the  large  ground  plane  simulat- 
ing the  infinite  natural  ground  must  have  a conductivity  of  40  a . Also, 
the  antenna  and  the  ground  plane  for  the  scaled  model  should  have  been 
fabricated  with  a material  whose  conductivity  was  40  times  higher  than 
the  copper  and  aluminum  used  for  the  monopole  and  ground  plane, 
respectively. 

For  practical  reasons,  copper  was  used  for  both  the  antenna 
and  for  the  finite  ground  plane.  Copper  mesh  screen  was  used  for  simu- 
lating the  infinite  natural  ground  of  the  reduced  model.  The  copper 
mesh  screen  used  was  in  the  form  of  a rectangular  screen  5 feet  by 
5 feet  in  size.  This  represents  6 wavelengths  on  one  side  at  the  upper 
end  of  the  frequency  band  and  1 wavelength  long  for  each  side  at  the 
lower  end  of  the  band. 

Due  to  the  finite  size  of  this  simulated  infinite- ground- 


plane,  the  radiation  patterns,  in  general,  tend  to  have  their  maximum 
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radiation  intensity  off  the  horizontal  axis,  where  ordinary  polar  angle 

A. 

(9)  in  the  spherical  coordinate  system  is  90  degrees. 

A quarter  wavelength  monopole  on  an  infinite  conducting 
ground  plane  has  a radiation  pattern  in  the  form  of  sin0 , which  gives 
the  maximum  intensity  of  the  radiation  field  in  the  direction  of  9 = 90°. 
Leitner  and  Spence  (Ref.  3)  have  showed, theoretically,  the  dependence 
of  this  pattern  on  ground  plane  size.  For  a circular  ground  plane  of 
radius  a,  the  far-field  radiation  patterns  were  plotted  for  ka  = 3,  4,  5 
and  6.  5.  They  all  have  maximum  intensity  off  the  horizontal  axis 
where  0 = 90°. 

Therefore,  the  results  obtained  from  this  experiment  are 
expected  to  differ  somewhat  from  the  theoretical  results  where  the 
ground  plane  was  assumed  to  be  infinite 

4.  3 Experimental  Technique 

The  radiation  patterns  obtained  in  this  experiment  are  taken 
in  an  x-z  plane  as  a function  of  polar  angle  9 , as  shown  in  Fig.  4.  1 . 
When  the  antenna  and  the  ground  plane  are  oriented  in  this  way,  the 
axis  of  the  monopole  coincides  with  the  z axis  of  a rectangular  coordi- 
nate system,  the  antenna  is  vertically  polarized.  This  pattern  is  some- 
times called  the  "E  plane"  pattern  where  E^  , with  0 constant  as  a 
function  of  9 , is  the  quantity  actually  measured. 

The  antenna  range  must  satisfy  a number  of  important  con- 
ditions before  it  can  be  used  for  antenna  radiation  pattern  measurements. 
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Fig.  4. 1.  Radiation  patterns  in  the  vertical  plane 

The  distance  between  the  transmitting  antenna  and  the  test  antenna  must 
be  sufficiently  large  to  assure  an  accurate  far-field  radiation  pattern. 
For  an  accurate  far-field  measurement,  the  antenna  under  test  must 
be  illuminated  with  a plane  wave  front.  Since  perfect  plane  wave  is 
only  possible  at  an  infinite  distance,  some  practical  limits  must  be 
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established  as  a guide  to  an  acceptable  distance  between  the  transmitting 
and  receiving  test  antenna.  A common  criteria  is  that  the  phase  differ- 
ence between  the  center  and  the  edge  of  the  antenna  under  test  be  no 
greater  than  A/ 16  (Ref.  11).  Then  the  minimum  acceptable  distance 
between  the  antennas  is  given  by  the  equation 

4a2 

R > ~ (4.1) 

where 

6 « R and  6 « a . (4. 2) 


The  distance  requirement  is  shown  in  Fig.  4.2. 


R2  + a2  = (R  + 6)2 


(4.3) 


6 


(4.4) 


The  range  used  for  this  experiment  was  constructed  with 
R ~ 20  meters.  Since  the  radius  of  the  finite  ground  plane,  a,  is  as 
large  as  A/4 


4a2 


A 


(4.  5) 


where  A varies  from  1.  5 meters  to  0.25  meter  from  200  MHz  to  1200 
MHz,  the  distance  R is  more  than  enough  to  satisfy  the  condition  (4.  1). 

Another  possible  source  of  error  in  the  radiation  pattern 
measurement  is  due  to  ground  reflection,  when  both  the  transmitting 
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Fig.  4.2.  Phase  difference  between  center  and  edge 
of  the  test  antenna 

and  receiving  antennas  are  located  close  to  the  ground  level.  One  way 
to  avoid  this  problem  is  to  utilize  a narrow  beam  transmitting  antenna. 
Another  possibility  is,  of  course,  locating  both  antennas  far  above  the 
ground  level. 

The  transmitting  and  receiving  antennas  for  this  experiment 
were  located  approximately  10  meters  above  the  roof  level  of  a three- 
story  building.  The  transmitting  antenna  used  was  a wide-band,  fre- 
quency independent,  log-periodic  sheet  triangular-tooth  structure 
designed  to  operate  above  150  MHz.  It  was  linearly  polarized  and  had 
a beam  width  of  approximately  60°. 

The  test  setup  used  for  the  experiment  is  shown  in  Fig.  4.  3. 
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Fig.  4.  3 Radiation  pattern  measurement  set-up,  block  diagram 
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4. 4 Measurement  Results 

The  '"  per ment  performed  here  was  designed  to  observe  the 
effect  of  a finite  ground  plane  upon  the  radiation  characteristics  of  a 
monopole  displaced  a certain  distance  above  an  infills'*'  conducting 
ground.  The  geometrical  arrangement  of  each  component  of  this  r.i- 
tenna  system  is  shown  in  Fig.  4.  4. 

The  values  used  for  a,  h,  and  £ in  the  experiments  were  as 
follows.  The  physical  length  of  h for  Fig.  4.  4 was  fixed  so  that  it  is 
one  quarter  wavelength  at  the  upper  end  of  the  test  frequencies,  which 
ranged  from  200  MHz  through  1200  MHz  for  a flat  disc  ground  plane. 

The  spherical  ground  plane  shown  in  Fig.  4.  4(b)  also  had  a fixed  length 
of  antenna  which  was  one  quarter  wavelength  at  the  upper  end  of  the  fre- 
quency band  ranging  from  246  MHz  through  1475  MHz. 

The  frequencies  chosen  for  the  spherical  ground  plane  were 
based  mainly  on  the  available  size  of  copper  half  spheres.  A quarter 
wavelength  at  1475  MHz  is  2 inches  or  5.08  cm.  A set  of  half-spheres 
available  was  2 inches,  3 inches,  4 inches  and  5 inches  in  diameter, 
and  at  the  frequency  of  1475  MHz,  they  correspond  to  A/4,  3A/8,  A/2, 
and  5A/8  in  diameter.  The  radiation  patterns  are  later  studied  as  a 
function  of  ground  plane  size  in  terms  of  wavelength  A . 

The  experimental  results  are  arranged  first  in  Fig.  4.  6 to 
see  the  effect  of  a ground  plane  size  when  it  is  located  at  a given  height 
(a)  above  the  infinite  ground  plane.  The  ground  plane  diameters  were 
changed  from  0.0156  meter  to  0.  125  meter  which  simulates  the 


gg  UH-J-s  TgRgPg  sew  WgV.Ji.igj1 
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Fig.  4.  4 Geometrical  arrangement  of  an  antenna  and  ground  plane 
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diameters  of  0. 625  meter  to  5 meters  of  the  original  antenna  system. 

These  diameters  correspond  to  X/ 16  through  X/  2 at  1200  MHz  and 
30  MHz.  The  changes  in  the  diameter  were  made  in  four  steps,  with 
the  size  of  the  ground  plane  doubling  its  diameter  each  time.  Thus,  the 
diameters  of  the  flat  disc  used  for  the  experiment  were  X/ 16,  X/  8, 

X/  4 and  X/  2 . 

Second,  in  Fig.  4. 10,  the  location  (a)  of  the  ground  plane 
was  changed  for  the  given  size  of  a ground  plane  diameter.  The  steps 
taken  for  different  a v’ere  the  same  as  discussed  above  for  the  variable 
si zt  o';  ground  plane,  namely,  a=X/16,  X/8,  X/4,  and  X/ 2 . 

In  Fig.  4.  5,  a monopole  above  a semi- spherical  ground 
plane  of  various  size  has  been  tested  for  its  radiation  patterns.  When 
the  ground  plane  size  becomes  larger  with  respect  to  wavelengths,  side 
lobes  become  noticeable.  When  ka  « 2 , the  side  lobe  level  shows  over 
60  percent  of  the  .najor  lobe  level.  However,  for  smaller  values  of  ka  , 
side  lobe  effect  is  negligible.  Typical  experimentally  observed  beam 
widths  measured  between  3db  points  range  from  approximately  30  degrees 
to  70  degrees.  These  patterns  are  later  compared  with  the  theoretical 
results  obtained  in  Chapter  5. 

Figures  4.  6 through  4.  13  .show  that  when  displacement  be- 
tween the  disc  ground  plane  and  a simulated  infinite  ground  plane  is  less 
thar  X/2  , major  lobes  are  directed  considerably  off  from  the  horizontal 
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position  with  a small  side  lobe  along  horizontal  direction.  This  type 
of  phenomena  can  be  also  observed  for  a center  fed  dipole  of  1. 5A.  in 
length  (Ref.  14).  This  indicates  that  by  using  a monopole  with  a finite 
ground  plane  displaced  a certain  distance  away  from  an  infinite  ground 
plane,  the  system  behaves  like  a longer  dipole  antenna.  A finite  ground 
plane  and  a gap  between  the  finite  and  infinite  ground  plane  simulates  an 
extension  of  an  antenna  length.  It  is  also  noticed  that  the  radiation  pat- 
terns generally  are  more  sensitive  to  the  gap  distance  a than  the  ground 
plane  diameter  D.  In  particular,  when  the  gap  distance  a = A./4, 
the  pattern  shows  a considerable  difference  fr  )m  other  patterns. 

In  summing  up  the  results  of  radiation  pattern  studies,  it 
can  be  concluded  that  the  finite  ground  plane  used  in  conjunction  with 
an  infinite  ground  below  causes  radiation  patterns  mainly  to  be  directed 
off  from  the  horizontal  direction  due  to  an  equivalent  antenna  arm  created 
by  the  finite  ground  plane  and  the  distance  to  the  infinite  ground.  Although 
minor  lobes  are  present  in  most  cases,  their  relative  levels,  compared 
with  major  lobes,  are  insignificant.  Also,  it  shows  that  the  distance 
a is  more  significant  for  the  radiation  patterns  than  the  diameter  of 
a finite  ground  plane  D . Finally,  there  seems  to  be  sc  ne  similarities 
in  radiation  patterns  between  the  monopole  with  hemispherical  ground 
plane  and  the  sleeve  dipole  (Ref.  15).  However,  a detailed  comparison 
is  not  possible  because  very  few  corresponding  dimensions  between  two 


antennas  are  available. 


Radiation  patterns  I Eg  I for  a monopole  with  a hemispherical  ground  plane  of 
various  diameters 


D*0625m 
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Fig.  4.  8.  Radiation  patterns  !E  I2  for  various  ground  plane 

0 

diameters  (D)  with  a = 2.  5m  , 
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Fig,  4,10.  Radiation  patterns  |E^f^  for  various  ground 
plane  locations  (a)  with  D = 0.  625m 


Fig.  4. 13.  Radiation  patterns  lE^i*  for  various  ground 
plane  locations  (a)  with  D = 5.  Om 


CHAPTER  V 


THEORETICAL  ANALYSIS* 

5. 1 Introduction 

In  this  chapter,  the  theoretical  model  analyzed  is  a 
monopoly  antenna  on  a hemispherical  ground  plane  that  is,  in  turn, 
above  an  infinitely  large  conducting  ground.  Far  zone  electromag- 
netic fields  are  computed  as  a function  of  hemisphere  radius  and  an- 
tenna length.  The  model  chosen  for  the  analysis  is  shown  in  the  Fig.  5. 1. 


Fig.  5.  1.  Theoretical  model  of  an  antenna  system 
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Using  image  technique,  the  original  model  is  divided  into 
two  parts  as  shown  in  the  Fig.  5. 1(b)  and  (c).  Each  configuration  is 
separately  studied  and  later  the  far- zone  fields  due  to  both  antennas 
and  the  induced  currents  on  the  spherical  surfaces  are  superimposed 
to  obtain  the  total  fields. 

The  Section  5. 2 deals  with  a development  of  a general  far- 
zone  fields  expression  from  a given  source  configuration,  and  the  Sec- 
tion 5.  3 includes  a far- zone  field  calculation  due  to  two  monopoles 
separated  by  the  diameter  of  the  sphere.  The  Section  5.  4 computes 
the  induced  current  on  the  spherical  surface  due  to  the  monopole  and 
its  image.  And,  the  Section  5.  5 shows  the  far- zone  field  expressions 
due  to  the  induced  surface  currents  on  the  sphere,  and  the  total  far- 
zone  fields  due  to  the  antennas  and  the  surface  currents  on  the  hemi- 
sphere. Some  of  the  results  from  the  numerical  calculation  in  the 
forms  of  radiation  patterns  are  given.  Section  5.  6 contains  results  of 
radiation  resistance  theoretically  obtained  for  a monopole  with  a hemi- 
spherical ground  plane. 

5.2  General  Far-Zone  Field  Expressions 

5.2.  1 Classical  Formulation.  Let  us  postulate  a new 
form  of  Maxwell's  equations  for  harmonically  oscillating  fields: 
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V x E = -jwpAH  - J 
0 m 


V x H = J + jutfgE 


V * J = -j u>p 
m J m 


V • J = -j  cop 


(5.1) 


where 

Jm  is  a fictitious  "magnetic  current  density" 

and 


Pm  is  a fictitious  "magnetic  charge  density" 

From  the  above  equations,  the  following  equations  are 

derived. 

v ■ (e0E)  - p V • (p0H)  = pm  (5.2) 


Now  define  the  electric  and  magnetic  fields  as  sums  of  two  parts, 
i.  e. , 


E = E' + E" 
H = H'  + H" 


(5.3) 


where  E'  and  H’  correspond  to  the  fields  arising  from  the  actual 
current  J when  *J  = 0 and  thus  satisfying  the  equations 

Vx  E’  = -jcop0H'  V • (e 0Ef)  = p 

(5.4) 


= J + jcoe  ' 


^ • (fi0H’)  = 


V x H’ 


0 
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Similarly,  E"  and  H"  satisfy  the  equations 


VxE"  = Jm 


V * (£qE”)  = 0 


V x H"  = ja*nE" 


v ■ tojn  - 


For  the  set  of  equations  5.  4 , a vector  potential 
and  a scalar  potential  , can  be  defined  as 


and 


M0  H'  = VxA 


where 


-ja>/i0co*  = V ' A 


E'  + jwA  = - V 4> 


Then  the  following  two  differential  equations  are  obtained. 


V?  A + k2  A = -P0 J 


V2  4>  + k2  4>  = 


€0 


k ■ a'  Mo€o 


(5.  5) 


A 


(5.  6) 


(5.7) 


(5.8) 


(5.9) 


where 
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The  solutions  of  these  differential  equations  are 


A = jDCT  G<  ^7  ) J(Rf)dV* 


V' 


* = 7 JJjGi-?-)  pm  dV' 


ho  v,  R' 


(5. 10) 


where 


j e-ifc1R  - R'l 

41  IR-R’I 


and 


(5.11) 


R and  R1  are  the  radial  distance  from  the  origin  of  the 
coordinate  system  to  the  observation  points  and  the  source 
points,  respectively. 

Fields  E’  and  H'  can  be  written  as 


E’  = -jw  l A + i-  v(v  * A)  J 
k2 


(5.  12) 


and 


H'  = 


if- VxA 


(5.  13) 


Similarly,  for  the  set  of  equations  5.  5 , a new  vector 

potential  A and  a scalar  potential  4>  can  be  defined 
r m r m 


Ill 


€nE"  = -VxAm 
u m 


“jajji-C  = V • A 
J 0 0 m 


m 


The  resultant  differential  equations  are 


V2  A + k2  A = -cnJ 
m m Om 


and 


+ k2  $ 
m m 


m 


The  solutions  of  these  equations  are 


* j;  /^G( f)  V5'>dV’ 

The  fields  E"  and  H"  can  then  be  expressed 

E”  = - — V x A 
e0  m 

H"  = -jto  I A + — V(V  • A ) ] 

J 1 m ,2  nr  J 


(5. 14) 


(5.  15) 


(5. 16) 


(5.  17) 


(5.  18) 


(5.  19) 
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Therefore,  the  total  E and  H fields  due  to  both  electric 
and  magnetic  sources  can  be  written  in  terms  of  vector  potentials  as 


E = E'  + E"  = - jw [ A +—  V(  V - A)]  - V x A (5.20) 

k2  60 


H = H'  + H"  = 


i-  VxA-iU[Am  + iv(V  AJ] 


(5.21) 


In  the  region  where  J = Jm  = 0,  these  equations  are 
simplified  to 


E = 


V x Vx  A — - 
k2  € 


V x A 


0 


m 


(5.22) 


H 


V x V x A + ~ V x A 

k2  m 


(5.23) 
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Consider  the  coordinate  system  shown. ahove  in  the  Fig.  5. 2 for  far- 
zone  field  R > > R’  and  |R-  R'l  = R-  R’  cos  a = 1R!  - (R 

A /s,  /\ 

If  the  direction  of  R'  is  denoted  by  8,0  where  R , and  R'  are 
the  unit  vectors,  then 


cos  a = sin  9 sin  6 'cos  (0'  - 0)  + cos  6 cos  9 ' 
With  the  approximation  shown  above 

IR  - R'l  » R - R'  cos  a 


and 


A = 


M0e 


-jkR 


4ttR 


ff  J'H R>ikR’  C°S  “dV' 


Similarly 


€oe 


■jkR 


m 


47tR 


J f j Jm(R')ejkR  c0SQ'dV' 


Then 


e~jkR 

5 ■ s <e>  »> 


e n e"jkR 

= -Q~\r  “ *> 


• R)- 


(5.24) 


(5.25) 


(5.26) 


(5.27) 


(5.28) 
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in  spherical  coordinate  system,  if  R > > A , then  all  terms  except 
jkN„  and  jkN . can  be  neglected.  Therefore, 

t7  0 


VxA 


Since 


V 


-jkR 


4irR 


0 


jk  [n  J - N„0 


1 ■ * 


V 


-jkR 


4irR 


/\  _ 
(Rx  N) 


(5.  32) 


VxA 


V x 


V 


-jkR 


4?rR 


N (0,  4>)  = Vx$N 


= -jk*  (R  x N) 


where 


* 


V 


-jkR 


4ffR 


VxVxA  = Vx(VxA) 

where  we  can  take  V x A as  a new  vector  ( W)  equivalent  to  N above. 


VxVxA 


V x W = (-jk)  (ft  x W) 


= k2 


/\  /\ 

(R  x (R  x N)  ) 


. V 


-jkR 


47rR 


(N  - NrR) 


(5.33) 
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where 


N-NrR  . Nt 


Finally,  the  far- zone  field  for  any  antenna  is  given  as 


E = jk  77q  [ -Nt  +1$  x Mt  J 


rJkR 

4ttR 


_ zs  _ P-ikR  , /v 
H = -jk[Mt  + RxNtJ  — j— = i®xE 


(5.34) 


4ttR  rj 


o 


These  expressions  show  that  only  the  transverse  components 
of  M and  N enter  into  E and  H far- zone  expressions. 

5.  2.  2 Stratton- Chu  Integral  Formulation . In  the  previous 
section,  a general  far-field  expression  of  E and  H have  been  developed 
from  Maxwell's  equations  with  a given  current  or  source  distribution. 

In  this  section,  the  same  far-field  expressions  of  electromagnetic  fields 
are  developed,  using  Stratton-Chu  integral  formulations. 

Let  P and  (5  be  two  vector  functions  of  position  with 
the  proper  continuity,  then 

J f fy  (Q’VxVxP  - P • Vx  V xQ)  dV 

(5.  35) 

= / / (PxVxQ-QxVxP)  • dS 
S 

where  S is  a regular  surface  bounding  the  volume  V . The  above 
integral  equation  is  the  generalized  Stoke's  theorem. 
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If  the  field  vectors  E and  H are  assumed  to  be  e ^ 
time  dependent,  the  Maxwell's  equations  can  be  written  as 


V x E - jtopH  = -7 
J ^ m 


V x H + ju>e  E = J 


V • H 


i 

— p 
p m 


V • E 


The  medium  inside  volume  V is  considered  to  be  homo- 
geneous and  isotropic.  The  quantities  J"m  and  pm  are  the  fictitious 
densities  of  magnetic  current  and  magnetic  charge. 


Z 


Fig.  5.  3 Notations  for  Stoke 's  theorem 
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We  consider  a volume  V bounded  by  the  closed  surface 

Stf  Sr  ®2  “ ' 38  s^own  *n  Fig-  3.  The  notation  S in  Eq.  5. 35 

means  a sum  of  SM  S+  . . . S . 

0*  I n 

From  Maxwell’s  equation,  expressed  in  Eqs.  5. 36,  the  pair 
of  vector  Helmholtz  equations  can  be  derived 


VxVxE-k2!  = jwpJ  - V x J m 


7xVxH*kJH  = jweJm  + V x J 


(5.  37) 


Where  currents  and  charges  of  each  are  related  by  the  continuity 
equations 


V ■ J - jwp  = 0 


V • J 

m 


0 


(5.38) 
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and 


k2  = o>2/ie 


In  Eq.  5.35,  let  P = 2 and  Q = 0$  where  £ is  an 
arbitrary  unit  vector  and  0 = e ft.  Distance  r is  measured  from 
the  element  at  (xf,  y’,  z’)  to  the  point  of  observation  (x,  y,  z);  i.  e. , 

r = ^<x-xf  + (y  - y')2  + (z  - z’)s  = IH-R'I 


(5.  39) 


R = (x2  + y2  + z2  ) 2 and  R' 


+ y 


>2 


2 


Then, 
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_ /\ 
V x Q = V0xa 


VxVxQ  = ak?0  + v(a  • V0) 


(5.  40) 


V/VxP  = k2  E + ju>  p J - V x J 


m 


Substituting  Eqs.  5. 40  into  5. 35  gives 

Iff  l 0a  * (k2  E + jwpJ  - V x J ) - E • (ak20  + V(a“*  V0)  ] dV 
y m 


= Jf  [ E x V 0 xli  - 0^a  x V x E ] * dS 
S 


(5.  41) 


Since 


E • V(sN  V0) 


V • U • V0)E  - (‘ a • V0)V  • E 


(5.  42) 


and 

VE  = 

Iff  v • U * V0)  EdV  = / f (a  • V0) E • dS 
V s 

= z • Jf  ( n"*  E)  V0  dS  (5.43) 
S 

Equation  5.  41  becomes 

Iff  [jwjtxJ 0 - Vx^m0  + | P 1 dV 
V 

= ff  [jwp(nxR)0  f(nxE)x  V0  + (n*  E)  V0  - .ixJ  0]dS 
S 
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where 


dS  = n dS 


Notice  that 

If  ( E x V0  x a)  • n dS  = //  n x (E)  x V0)  • a dS 
and 

//(0a  x V x E)  • ndS  = / / (0 a x (jo>pH  - Jffl)  • n dS 

= *//jw/i(Hxn)0  • adS-  //  nx  Jm0  . adS 

Equation  5. 41  can  be  further  reduced  to 

Iff  (jun J0  - J X V0  + i P V0)  dv 
y mg' 

= //  f jcopt(n  x H)  + (nxE)  x V0  + (n  . E)  V0  ] dS 
S 

(5.  44) 

with  the  help  of  identities: 


///  V x (Jm0)  dV  = Iff  V0  x Jm  dv  + ///  0 V x Jm  dV 


and 


(5. 45) 


///  v x (Jm0)  dV  = //  n x Jm0  dS 
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Thus 


-/J7*vx  Jn/V  = -//'nx  Jm0dS- ///Jmx  V0dV  (5.46) 

In  Fig.  5. 3,  the  surface  S includes  Sj,  Sg,  ...  and 
a volume  surface  surrounding  the  observation  point  (x,  y,  z). 

Since 

eJkr  j /\ 

V0  = V(^— ) = (4-  - jk)-£-j-?0  (5.47) 

where  r^  is  a radial  unit  vector  pointed  toward  the  center 
of  the  spherical  surface  . Thus,  r^  = ft  for  the  surface  . 

Evaluating  the  surface  integral  of  the  Eq.  5. 44  over  the  surface 
Sq  , notice  that 

(n'x  E)  x'n't  (n>*  E)  n = -/n'x  (n'x  E)  + (rN  E)  n 

= -'n(E  .li)  + E(n'*/n)  + (n  • E)^  (5.  48) 

= E 


Also,  when  the  radius  of  the  sphere  r^  represented  by  Sq 
is  reduced  to  zero,  the  only  terms  remaining  from  the  surface  integral 
are  terms  involving  l/r^2  : 
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IJ  £ (n^x  E)  x V0  + ju)fjL(nx  H)p  + (n'*  E)  V0  J dS 
S0 


2*  r 

I I 

0 0 


('r'.x  E)  x'r,  — — — 


v i*  ' «« 


0 


0 2 

ro 


r^2  sin  6 dd  d0 


where 


*(tt) 


represents  other  terms  involving 


As  Tq  approaches  zero,  the  integral  uecor.ies 

J J [(n'x  E)  x V0  + ju)p(n  xi)0  + <n-  E)  V0]  dS’ 

so 


4ffE(x,y,z) 


Equation  5.  50  then  becomes 


E(x,  y,  z)  = J J J (jwp  J0  - J m x V0  + — p V 0)  dV' 


-J-  /'/ 

4tt  j • 


[ jcup('n'x  H)0  x (n  x E)  x V0  h 


If  the  volume  does  not  contain  any  sources,  Eq.  5.  5 


■) 


(5.  49) 


(5.50) 


(I:-  E)  V0]dS’ 
(5.51) 

further 


reduces  to 
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E(x,  y,  z)  = - -jf 

[ju>/i(frx  H)0  + (frxE)  xV0+(fr-  E)  V0]  dS*  (5.52) 

Also,  if  all  currents  and  charges  can  be  enclosed  within  a sphere  of 
finite  radius,  then  the  field  is  regular  at  infinity  and  either  side  of  S 
may  be  chosen  as  its  interior.  The  surface  integral  of  Eq.  5.  51  rep- 
resents the  contribution  of  sources  located  outside  S . If,  therefore, 

S recedes  to  infinity,  the  contribution  from  those  sources  vanishes. 
Discarding  the  fictitious  magnetic  charges  7m  , Eq.  5.  51  becomes 

E(x,  y,  z)  ^ -TR-Jff  f + V0]  dV'  (5.53) 

Equation  5.  52  is  the  field  that  would  be  produced  from  surface  S by 
electric  current  density  K,  by  magnetic  current  density  Km  , and  by 
surface  electric  charge  of  density  or  ; where 

S 

ftxir  - -K,  'n'x  E = K or  = -e'fr-  E (5.54) 

III  b 

If  Eqs.  5.  52  and  5.  53  are  applied  to  the  far- zone  where 
0 - e^r/r  , then 

jkR  - jk&  • E’ 

, _ e 
0 = 


R 


(5.55) 
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and  Eqs.  5. 52  and  5. 53  can  be  written  as 


E(x,y,z)  = iff 


•t-  jj  [ja>ji(nxH)0  + (nxE)xV0]  dS' 

CO  c ^ 

Sn 


(5  56) 


and  if  Sp  Sg  . . . Sn  are  perfectly  conducting  boundary,  then  n xE  = 0 


Therefore, 


E(x,y,z)  = 4?  ff 


SpS2,  ••  • sn 


[jw/i(n  x H)0  ]t  dS’  (5.57) 


where  subscript  t signifies  transverse  component  to  the  radial  direc- 
tion. The  proof  of  these  equations  5.  56  and  5.  57  is  given  in  Appendix 
B.  Of  course,  Eq.  5. 56  is  equivalent  to  Eq.  5.  34  in  Section  5.  2. 1. 
Similarly,  it  can  also  be  shown  that  in  the  far- zone, 


H(x,  y,  z)  = ^ ff 


Sp  Sg.  • • , S 


[jo>e(n  x E)0  +(nx  H)x  Vd^dS' 


n 


(5.58) 


5.  3 Far- Zone  Field  Expressions  for  Two  Linear  Antennas 

Consider  an  antenna  and  its  image  with  a sinusoidal  current 
distribution  as  shown  in  Fig.  5.  4. 
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N = R')ejkR,’Rdv’ 


when 


— /\ 

R * R = R*  cos  a = R'  cos  6 cos  Q ' + 


R - Unit  vector  in  radial  direction 


sin  9 cos  ( 0 - 0’) 


In  the  problem  considered  only  line  sources 


are  present.  Therefore, 


N = / I (R')ejkR 


where  R’  = z*  and  R’  * R = z’ cos  0 for  antenna  A 
and  R - - z z and  R ’ • & = +z>  cos  q for  antenna  B 


Nz  " / Imsink(d-z-)ejk2'cosedz,  + 7I 

P J r 


sin  k (d  + z’)  e^kz'  cos  Gc 


Changing  the  variable  for  the  second  integral  by  z’  = -z” 
the  integral  becomes 


(5.61) 


-i  ” Sm  k(d  ’ Z<>  ^ C°S  9 dz’  = / I sin  k(d  - z")  e'jkz"  c°s  6 d,. 


(5.  62) 


= <„,/  smk(d-z’)!eikz' cos9  + e'ikz’cos9|  dz' 

Cl  J 


= 2 *m  .1  sin  k(d  - zf)  cos  (kz*  cos  9)  dz’ 


(5.  63) 
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Since 


i [ sin  (a  + b)  + sin  (a  - b)]  = sin  a cos  b 


N = 


d 

I / sin  [ k(d  - z')  + kz’  cos  0]  + sin  [ k (d  - z’)  - kz 
a 


’ cos  0]  | dz’i 


= Im  / | sin  [ kd  + kz'(fc0s-l )]  + sin  [ kd  - kz‘(l  + cos  0)]  j dz’ 


(5. 64) 


Also,  using 


r 1 

J sin  (a  - bx)dx  = -r—  cos(a  - bx) 


Nz  = lm  [k(L  - cos-9)  cos[  kd  - kz'O  - cos  9)]  + k(l  VSbTSJ 


cos  [kd  - kz'(l  + cos  6 )]  J 
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N = I COS  (kd  COS  9)  [ -rrr — -rr  + Y-K - jrr  1 

z m ’ 1 k(l  - cos  9)  k(l  + cos  9)  1 

r (cos[  k(d  - a)  + ka  cos  91  cosfk(d  - a)  - ka  cos  9] ) 
mj  k(l  - cos  9)  L + k(l  + cos  9)  L j 


2 I cos  (kd  cos  9)  I 
m m 


k sin2  9 


k sin2  9 


£(1  + cos  9)  cos[  k(d-a)  + ka  cos  9] 

+ (1  - cos  9)  cos  [ k(d-  a)  - ka  cos  9j 


Using  trigonometric  identities,  Eq.  5.  65  can  be  written  as 


(5.  65) 


2 I cos  (kd  cos  9)  21 

N„  = r — — [ cos  k(d  - a)  cos  (ka  cos  9) 


k sin  9 


k sin  9 


- cos  9 sin  k(d  - a)  sin  (ka  cos  9)] 


2 I 


m 


[ cos(kd  cos  9)  - cos  k(d  - a)  cos(ka  cos  9)  + cos  9 sin  k(d  - a) 


k sin2  9 

sin(ka  cos  9)] 
(5.  66) 

When  a = 0,  the  expression  reduces  to  the  radiation  vector  of  a dipole 
which  is 


2 I 


N = 
z 


m 


k sin  9 


[ cos(kd  cos  9)  - cos  kd] 


(5.  67) 


k 
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The  6 component  of  N is 

•2  I 

Ng  = -Nz  sin  0 = -g  [ cos(kd  cos  0)  - cos  k(d  - a)  cos(ka  cos  0) 

+ cos  0 sin  k(d  - a)  sin(ka  cos  0)] 


The  far-zone  electric  field  is,  from  Eq. 


E = jkr70  [-Nt] 


5.34, 


e-ikR 

4ttR 


(5.68) 


(5.34) 


when 


Mt  = 0 


Therefore,  the  electric  field  due  to  two  monopoles,  E , , becomes 


E0  = 


_ -jkR  p 

J ^ome  cos(kd  cos  0)  - cos  k(d  - a)  cos  (ka  cos  0) 


2ttR 


L 


sin  0 


| cos  0 sin  k (d  - a)  sin  (ka  cos  9) 

sin  0 J 

(5.  69) 
and 


H = — § x E = — En# 

% Q 


(5.70) 
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The  radiation  patterns  due  to  the  two  monopoles  discussed 
above  are  the  plots  of  angular  dependent  terms  of  E as  a function  of 

e . 


•p(Q)  - c°s(kd  cos  9)  - cos  k(d  - a)  cos(ka  cos  6)  + cos  0 sin  k(d  - a)  sin(ka  cos  0) 


sin  0 


(5.71) 


When 


a = d 

p(9)=  cos(kd  cos  6)  - cos  kd  (5.72) 

sin  8 

This  expression  is  the  angular  term  of  a far-zone  field  expression  of  a 
dipole. 

5.  4 Induced  Current  on  a Spherical  Surface 

5.  4.  1 Induced  Current  on  a Conducting  Sphere  Excited  by 
a Monopole.  The  induced  current  on  a conducting  sphere  excited  by  a 
monopole  antenna  erected  on  its  surface  was  treated  by  Papas  and  King 
(Ref.  16).  The  current  along  the  antenna  was  assumed  to  have  a form  of 
sinusoidal  distribution  and  was  assumed  to  be  independent  of  the  current 
on  the  sphere. 

This  solution  is  reviewed  here,  because  the  same  method 
will  be  adopted  to  calculate  the  induced  current  on  the  spherical  surface 
due  to  the  image  antenna  located  along  the  negative  Z axis. 


,^i js.*  ■ v .u.  j,  „«  j ^ v v i"-'  v, . i >■.■  ■*  vrv&i&!*xzFrvf*rrTr  i w* : mv^u't*1  y<  ~mj.WfV!.T-' 
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Z 


Fig.  5.  5.  A monopole  antenna  above  a spherical  ground 

plane. 

From  the  general  reciprocity  theorem  (Appendix  A) 

J J [ 'n • (Ej  x H2  - E2  x Hj)  ] dS  = 0 (5.  73) 

s 

where  and  Hj  are  the  electromagnetic  fields  associated  with  a 
current  density  J ^ and  E g,  Hg  with  J g.  Also  s is  the  closed  surface 
of  the  volume  of  empty  space  bounded  by  the  surface  of  the  antenna  and 
the  sphere  and  an  imaginary  boundary  at  infinity.  'n'  is  an  external  normal 
unit  vector  of  this  surface. 

The  assumption  made  here  is  that  the  two  sets  of  fields, 

Ep  Hp  and  Eg,  Hg,  are  defined  in  the  free  space  and  are  bounded 
by  the  same  geometrical  surfaces,  but  not  necessarily  satisfying  the 
same  physical  properties. 
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Let  fields  E ^ and  be  caused  by  the  actual  currents 
in  the  antenna  located  at  d = 0 and  a R < a + h as  shown  in 
Fig.  5.4  and  the  current  on  the  surface  of  the  sphere. 

The  boundary  conditions  for  fields  and  are: 

n x Ej  = 0 on  the  spherical  surface  Sg 

ft  x Hj  = Klg  " M ” M Sg 

The  current  on  antennas  I ^ and  K^g  on  the  sphere  is 
assumed  to  be  maintained  by  an  appropriate  generator  or  a distri- 
bution of  generators  in  the  antenna  in  such  a way  that  no  current 
exists  in  the  0 direction  of  the  antenna  or  on  the  sphere. 

Therefore, 

Kls(0)  = focis(0)  (5-74) 


and 


2zr 

I ve>  a sin  0 d0  = ljg(0) 


(5.75) 


Ilg(0 ) is  the  total  current  crossing  a parallel  of  latitude  on 

the  sphere. 

Also,  the  total  current  in  the  antenna  at  6 = 0 is  in  radial 
direction  in  the  spherical  coordinate  system.  Due  to  the  rotational 
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symmetry  about  the  axis  .9  =0,  the  line  integral  of  the  tangential 
component  of  around  the  surface  of  the  antenna  is  equal  to  the 
total  axial  current  I^(R)  in  the  antenna. 

Thus, 


2tr 


6 H.  - dl  = / H 
1 0 


10  R0  d0  “ 


/ / 


(cross  section  \ 
of  the  antenna J 


llA 


dS  = 


■ iia(r) 


(5.76) 


The  second  set  of  fields  Eg  and  Hg  are  not  the  actual 
ields  that  are  related  to  the  real  currents  i ^ in  the  antenna  or  on 
the  sphere.  They  are,  instead,  fictitious  fields  that  satisfy  the  same 
field  equations  as  the  first  set  and  that  must  be  defined  over  the 
same  geometrical  surfaces.  Since  fields  Eg  and  Hg  are  not 
required  to  be  defined  on  the  boundaries  with  the  same  physical 
properties  asset  Ej  and  Hp  let  the  volume  possessed  by  the 
antenna  for  Ej  anc  Hj  be  empty  space  for  Eg  and  Hg.  Also, 
on  the  surface  of  the  sphere,  Sg,  fields  Eg  and  Hg  are  required 
to  possess  the  following  prescribed  form  by  some  appropriate 
set  of  generators: 


E 

H 


20 

26 


0 

0 


U sin  6 P (cos  9) 
n 

(5.77) 


0 


The  entire  field  in  empty  space,  therefore,  has  the 
components  EgR,  Eg^  and  H,,0  which  must  satisfy  the  field 
equations  and  the  prescribed  boundary  conditions. 


i irrftiHfifcti 


.fin  rv-.T tier  r iTU 
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Returning  to  Eq.  5. 73, 


///[  n * (Ej  x H2  - E2  x Hx)  ] dS  = 0 (5.73) 


the  surface  S can  be  divided  into  two  parts  S ; the  surface  of  the 

s 

sphere,  and  SA,  the  surface  of  the  antenna. 


JJ  [ n • (Ej x H2  - E2  x Hj)]  d6  = ff[  (nx  Ej)  + E2  • (n x Hj)]  d8 


+ jQT  [ H2  • (n  x E^  - Hj  • (nx  E2)j  dS 
SA 

(5.78) 

Now  since  n * 4 over  SA  and  since  n = -R  over  S and 

s 

with  the  boundary  conditions  described  above,  each  integral  can  be 
rewritten 

a + h 2n 

//{h2  • (S'*  Ej)  - Hj  (px  E2)}dS  = EmdR^  H20  r0  d0 


a+h  2n 

'L  ***  • l Ro 110 

(5. 79) 

The  volume  occupied  by  the  antenna  is  empty  space  for 
H20  • Therefore,  for  a thin  antenna  where  the  radius  RQ  goes  to  zero, 
the  line  integral 
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Therefore, 


£ V**  - 0 


a+h 


(8  x Ej)  - Hj  (8  X E2)|d8  = - / E2R  Iu(R)  dR 

(5.80) 


The  surface  integral  on  the  sphere  becomes 


III  H2  • (8  x Ej)  + Eg  • (8x  Hj)]  dfi  = //  E2  • dS 

(5.  81) 

because  of  the  ooundary  condition  n x E ^ = 0 on  Sg. 

From  Eqs.  &.  75  and  5. 77 

II  V = U29  yejadO 

s o 

s 

Therefore, 

JJl  n • (E1  x H2  - E2  x Hj)]  dS  = 0 (5.73) 

s 

which  implies  that 


a+h 


/E29Ils(9)ad«  ‘ / 


E 


a 


2R 


I^G?)  dR 


(5.  82) 


This  is  the  integral  equation  to  determine  I^g(0 ) knowing 
IU(R)  and  the  prescribed  field  E^  will  assume  that  the  surface 
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current  I^g(0 ) takes  the  form  of 

00 

lls(0)  = Z BnPn(cos6')  (5.83) 
n=0 

where  the  coefficients  are  to  be  evaluated  from  the  integral 
equation. 

For  the  second  set  of  fields  Eg  and  Hg,  the  volume 
surrounded  by  the  closed  surface  S , S.  and  the  boundary  at  infinity  is 
empty.  Furthermore,  Eg^,  Egfl  and  Hg^  are  the  only  non-zero 
components  due  to  currents  on  the  sphere.  For  TM  fields  of  this  type, 
the  fields  are  found  from  the  scalar  wave  equation 

V2(-g-)  + k!  (-g-)  = ° (5.84) 


where 


k * “VVo 

The  solution  of  the  equation  is 

oo 

“ = L AmPm(C0S  6 ’ Pm(kR)  <5-85> 

m =0 


where 


Pm(cos  6 ) is  the  Legendre  polynomial  of  order  m 
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Pm#®)  = W Hm Vl  <“*>  “ <“)  h®  0®) 

H^(kR)  : Hankel  function  of  the  2nd  kind 
m 

(2) 

h m(kR)  : Spherical  Hankel  function  of  the  2nd  kind 


The  components  of  the  fields  are  obtained  from  solution  u as 


E2R  “ [ k*  + ~T  ] u 


3R' 


F - i /JL \ 

E29  “ R 3R  30  U 


(5.  86) 


H = ft  ) ( gu  1 

20  * n M ' 


From  Eq.  5.  85  and  5.  86  , on  the  spherical  surface  R = a, 


E0fl  = -( ) LAP  (cos0)p  (ka) 

2d  a ' Lj  n m m ,rm 

m=u 


(5.  87) 


Also,  we  demand  that  E ^ take  on  the  spherical  surface: 


E29  = ^ sin  ® Pn(cos  9)  = U 


Pn+1l(cos  9)  - Pn_  jl(cos  9) 


2n  + 1 


(5.  88) 
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Equating  the  two  equations 


5.  87  and  5. 38.  leads  to 


aU 

* (2 n + 1)  ‘ k • Pn  “ ^(ka) 


a U 

(2n  b 1)  * k • pn  _ i (ka) 


(5.  89) 


Therefore, 


IV  + 02  1 

‘ a Pn  _ j(cos  0)  pn  _ 1(kR) 

aP  ,(cosQ)p  .(kR)  * 
n + r ' rn  + V 

L 3R>  J 

. (2n  + 1)  • k • pn  _ (ka) 

(2n  + 1)  * k ' pn  + L*(ka) 

(5.90) 

Assuming  that  the  antinna  current  takes  the  sinusoidal  form 


(R)  = I 


max 


sin  k(d  - R) 


(5.91) 


it  can  be  shown  that^Vppendix 


a + h 

c 

J 

a 


E2R  W*'  * 


U a I fp  - l(kd) 
max  n 

2n+  1 [pn  _ ^(ka) 


- cos  kh 


+ cos  kh 


Pn  . 

p . (kd) 

*n  + r 

Pn  _ y (ka) 

p ,'(ka) 

1 n + 1 

j , (ka) 

n + 1 ‘ 

i '(ka) 

n + 1 J 

(5.  92) 

along  0 = 0. 
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Also 


/ E29  'is'5 * * 8)  ad8 


It  00 

fE9Q  2 B P (cos  9)  ad© 
0 20  n=0  n 11 


B 

n 


U 2a 
2n  + 1 


(5. 93) 


Finally,  the  integral  equation  produces  the  result 


I 

max 

Vi(kd) 

Pn  . ^ka) 

P ,(kd) 

*n  + r 

2 

/n- 

c 8 kh  pn  -T  (ka) 

vTT® 

,/n*  l(ka) 

+ C0Skh^77Tki) 

(5.94) 

The  surface  current  Ijg(>)  is  now  determined  in  terms  of  the 

coefficients  B ; 

n 

00 

Ils(0)  = L BnPn(cos  9)  (5.83) 

n=0 


5.  4.  2 Induced  Current  Due  to  the  Image  Antenna.  In  the 

Section  5.  4. 1,  the  induced  current  on  the  spherical  surface  due  to  the 

sinusoidally  distributed  current  on  the  antenna  was  found  in  the  form  of  a 

modal  current.  The  coefficients  B of  the  infinite  series  that  gives 

n 

the  surface  current  were  evaluated  in  terms  of  known  parameters.  In 
this  section,  the  contribution  from  the  image  antenna  to  the  total  surface 
current  on  the  spherical  ground  plane  is  obtained  following  the  method 


WfWWJSW  IJ  ►'  JSW wn ^’7n?<^7Tr^? } ,*■*.'  -> ssjww’*  % j u?w  ■ w^w^yj^ST^^Rv^ii  -i,  5T WE?*^3WCS£j p^r^v  »p?q,  t ;vTr? v»t>  ^7T?>l*v?y7?gP7T^~^ft"*^  ^ vnr^?»jgy,r>’  7 v*  >.-w*'  c «? ■*'  j ^ 
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used  in  the  previous  section. 


Z 


Fig.  5.  6.  Image  antenna  with  a spherical  ground 

plane. 


The  geometrical  arrangement  of  the  image  antenna  and  the  spherical 
ground  plane  is  shown  in  Fig.  5.  6.  The  center  of  the  spnere  is 
located  at  the  origin  of  the  spherical  coordinate  system. 


L 
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The  electromagnetic  fields  and  the  currents  are  all  expressed 
in  the  same  way  as  they  were  in  the  previous  section  except  for  the  prime 
notation  to  distinguish  the  image  antenna  problem  from  the  real  one. 

From  the  reciprocity  theorem, 

N (E^  x Hg  - Eg  x Hj)  • dS  = 0 (5.95) 

where  surface  S is  now  the  closed  surface  enclosing  the  volume  bounded 
by  the  spherical  surface,  the  surface  of  the  image  antenna  and  a boundary 
at  infinity.  The  surface  normal  unit  vector  'n'  is  the  same  as  was  in 
Section  5.  4. 1 on  spherical  surface  Sg  . However,  on  the  surface  of 
the  image  antenna  the  unit  vector  is  in  the  positive  'O'  direction. 

The  integral  equation  5.  82  is  now  modified  to  become 

7r  , , ath  , , 

/ E26Ils(S)ad9  = ■ / E2R  riA(R)  ^ <5-  96) 

0 a 

Let  the  surface  current  Ijg  ( 0 ) on  the  sphere  due  to  the 
image  antenna  also  take  a modal  form  with  coefficients  of  the  infinite 

T 

series  B . Then, 
n * 

i 20  t 

Ils  (0)  = L Bnpn(cos  (5.97) 


Following  the  steps  shown  in  Section  5.  4. 1 


U[  k2 


aVl(cos  9>Pn-l(kR)  aWcos0W  + i<kR>, 
72n  + 1)  ■ k -pn_  i'tkaj  (Zn  + !)■  k #>n  + 1'(ka)  J 


(5.P8) 
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Since 


P .(cos  0)  = (-l)n”  1 at  0 = ir 
n-  1 


(5.  99) 


along  the  image  antenna 


Pn  + 1(eos9)  = (-l)n+1  at  9 = » 


(kR) 


E ' - ur  Ic2  i ^ I r t\n~  1 P.  n-V  /1\n  + 1 Pn  + 1 

E2R  +3rJ  j[a(-l)  ^;1)kpn  i,(ka)  -<-l) 


,(kR) 


k Pn  + l'<ka)! 


n-1  r 2 

r*  4 TTf  b-2  o if  n-1 


^ i(kR) 


= (-1)  Uf  k2  + _iL_  1 r n-i  n+  V 

1 ak2  J T2n+ 1)  kp  ^’(ka)  " aT^i  + l)  • k-  Pn  + 1’(ka)  J 


(5. 100) 


■/h(-DnE2iiI1A'dR  . (-1)" 
a n - lv 


- cos  k(d  - a) 


,'Wka)  P„  + !<“) 


Pn  + 1(ka> 


viM  'pTT7®)  + C0Sk(d-a)^krJ 


Pn+i'(ka)  (5.101) 


Since 


a+h  a+h 

J E2R' IlA,(R)dR  = '/  *29  V»)  ad9 
a R=a 


(5.  102) 
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where 


'll'W  - W sta  k<d  - R> 


and 


V(8>  ■ 1 Bn'  Pn(cos  #) 


(5.  97) 


Also 


U sin  0 Pn(cos  0)  = U [ 


Pn  + ^(cos  0)  - PR  _ (cos  0) 
___ 


(5. 103) 


Thus, 


n co  T _ 

J E ’ T B ’ P (cos  0)  • ad0  = B ’ -% 

o 29  i£o  n n n (2n  +1) 


(5.104) 


Using  the  orthogonality  relationship  for  the  Legendre  function 


/ Pm(  M)  pn(M)  “ i 2n  + 1 1 5m,  n 


2B„' 


■ /E'2RIlA'(R)dR  ' 


n + 1 


(5. 105) 


pn . i<ka> 


Pn  i(^)  »"n  _ p' — ' 

wi?7r^)-co£'i<d-a)v^5 
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where 


I p ,(kd) 
„ _ max  f - 1 

n " ~2  Lp— [W 


cos  k h 


Pn  . !<ka) 


Pn  + jOcd) 


On  - l‘(ka)  " Pn+  l'<ka) 


+ cos  k h 


pn+l(ka)  ] 


pn+  1 
(5.  110) 


Finally, 


1S^>-  W s 

n=odd 


l 


pn  - l(kd) 
pn- 


- coskh 


n-  1 


(ka) 


n + 1 


(kd) 


pn-  l'(ka)  ' pn+  l',ka) 


Pn  + l*ka^ 

+ coskh  p — rtirar]  pn(cos  e) 

n + 1 


(5.  Ill) 


5.  5 Far- Field  Expressions  Due  to  a Spherical  Surface  Current  Distribution 


The  radiation  vector  N due  to  an  electrical  source  I*  is 


found  to  be 


N = jfjfjf  I'ejk  R'  cos  ^ dv’  (5.  112) 

v 

where  \j/  is  the  angle  between  the  radial  lines  to  the  source  point  and 
the  observation  point. 

The  far- zone  electromagnetic  fields  are  obtained  from 


radiation  vector  N: 
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,i "H  e’ik  R 
4zR  6 


E.  =JSH-e-jkRN 


0 4irR 


0 


(5. 34) 


j qos  t \t 

Since  R'  = a at  the  surface  of  the  sphere,  the  term  e v in 

Eq.  5. 112  can  be  written  as 


jk  R’  cos  if/  _ jk  a [cos  0 cos  0’  + sin  9*  sin  0 cos  (0  - 0’)  ] 

(5. 113) 

The  current  density  I*  inside  the  volume  integral  is  the 
surface  current  per  unit  width  and  flowing  in  0"  direction.  Therefore, 

I'  = k(0')0S  = k(O')  [ cos  O’  cos  0'x1  + cos  0'  sin  0’  y' - sin  ©•'z] 

(5.  114) 


where  x,  "y,  and  z"  are  unit  vectors  of  the  Cartesian  coordinate  system. 
The  radiation  vector  N can  be  written 

_ 2n  ir  A 

N = f f k(0’)[  cos  0'  cos  0’  x + cos  0'  sin  0’  y - sin  0’ z] 
o o 

ejka[  cos  0'  cos  0'  + sin  0’  sin  0 cos  (0  - 0')]  ^ gin  Q,  dQt  d0t 

(5.  115) 

Breaking  N into  three  components 
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N x = J J k(®')  cos  O'  COS  0' 


ejka[  cos  0 cos  O'  + sin  0’  sin  0 cos  (0  - 0’)]  ^ gin  #f  dQf  d0t 


2v  i 

N = J J k(0')  cos  O'  sin  0' 
y oo 

ejka[  cos  0 cos  O'  + sin  O'  sin  0 cos(0  - 0’)]  ^ gin  #f  dQt  d0, 

(5.  U6) 

2v  ir 

N = -J  J k(0')  sin  O' 
o o 

ejka[  cos  0 cos  0'  + sin  0’  sin  0 cos  (0  - 0’)]  gin  Qf  dQf  Q0f 


Also 


2 if 

J k(0')  a sin  0’  d0'  = 1(0')  (5.  76) 

o 


Since  k(O')  is  independent  of  0' 


k(O')  = 


I(Q') 

2jra  sin  O' 


n=odd 

y Bn  Pn(cos  O’) 
jra  sin  O' 


(5. 117) 


Using  the  identity(Ref.  14,  p.  407) 


gjk  R cos  \f/ 


0)n(2n  + 1)  i 


(kR)  Pn(cos  4/) 


(5. 118) 
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where 


J„4<kR> 


(5. 119) 


and  expanding  Pn(cos  if/)  into  a finite  series  of  the  form 


c n 

Pn(cos  if/)  = -=-  P (cos  8)  + Yi  (cm  cos  + dm  sin  m0)  (c°s  8) 

m=l  1 ! 

(5. 120) 


and  evaluating  the  coefficients  c and  d^  by  using  the  orthogonal 
relationship  of  the  Legendre  polynomials  such  as 


2ir  ir 

/ f Pn(oos  *)  Pnm(cos  0)  cos  m0  sin  0 d0  d0  = (n-mj*  cm 

o o 

(5. 121) 

and 


27r  jt 

/ /Y  (0,0)  P (cos  0)  sin  8 d6  d0 
o o n n 


4tt 

2n  + 1 


[ Yn(8,  0)] 

0=0 


(5.  122) 


where 


Yn(8,0)  = Pn(cos  0)cos  m0  (5.  123) 


cos  if/  = sin  9'  sin  9 cos  (0  - 0’)  + cos  9 cos  0 


it  can  be  shown  that 
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2v  v 

I I Pn(cos  P^(cos  0)  cos  me*  sin  0 dO  d0 
o o 


2~y  P™(cos  O’)  cos  m0* 


(5. 124) 


The  coefficients  c and  d are  evaluated  to  be 
m m 


c 

m 


d 

m 


o(n  - m) ! 
^(n  + m)! 


o(n  - m)! 
4 (n  + m)! 


P^(cos  0')  cos  m0* 
P™(cos  0’)  sin  m0’ 


(5.  125) 


Therefore, 


Pn(cos  ip)  = Pn(cos  0')  Pn(cos  0) 


+ 2 Z fo  V Pnm(coS  0')  p”(oos  0)  cos  m(0  - o’) 
m=lv  ' ’ 


n 


(5.  126) 


2ir  it  -j  B2C  + 1 P2f+1  <cos  e,) 

«W  / 1 


o o 


2rr  a sin  0’ 


cos  0’  COS  0’ 


ikaf  cos  0 cos  0’ + sin  0 sin  0’ cos  (o  - o’)l  * „ 

e J a sin  9’  d9’  do’ 


Since 


(5.  127) 


ejka  cos  \p 


Tj  * 

Z 0)  (2n  + 1)  Y ■ (kal 
n=o  u 


P ( cos  0’)  P (cos  0) 
n n ' 


+ 2 


v (n  - m)i 
m= l{n  + t^ 


P™  (cos  0’)  P™  (cos  0)  cos  111  (o  - o’) 


(5.  128) 
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and 

cos  m (0  - 0')  = cos  m0r  cos  m0  + sin  m0  sin  m0* 

the  contribution  of  the  first  integral  with  respect  to  0* 
reduces  the  double  integral  to 


N*  ■/ 


i ? B2I+1  P2I+1(cos  0,) 


wa 


f )n  (2n  + 1)  jn(ka) 
n=o  u 


2-^— 4jf  P„  (cos  O')  (cos  0)  it  - cos 0 a2  cos  9’  d9’ 


^n(ka)  P^  (cos  O’)  P^  (cos  9)  cos  9’  d9’ j 


COS  0 


because 


2n 

f cos  m0’  cos  0'  do'  = 0 for  1 

o 

= for  m = 1 


(5. 129) 


J sin  mo’  cos  O'  do'  = 0 for  all 


m 


Similarly,  for  the  evaluation  of  , we  obtain 


N = 

y 


A 


Z b21+1  p2f+i(cos  c,)  ? ini|rnT 


n 


9 

2 ^(ka)  P^  (cos  9')  P^  (cos  9)  cos  9’  d9’ 


sin  O 


(5.  130) 
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If  we  ?et 


71 


= a 


77 

/2b 

o f 


2i+l 


P2£+1(C0S  Q,)  S 3n 

n 


2n  + 1 
n(n  + l) 


then 


# 

2 ^n^  Pn  ^C0S  Pii  ^cos  0)  cos  ®*  d0, 

(5. 131) 


Nx  = 7tC0S  0 
Ny  =7lsin0 


(5. 132) 


/\ 


Finally,  the  z component  of  the  radiation  vector  N can 


be  written  as 

2tt  2ti  r 

N = - / / k(0  ')  sin  0 * ejka  l cos  6 cos  6 ' + sin  e sin  0 ' cos  0?  ]a  sin0  ’d0  'd0r 
Z 0 0 

0 0 


7ra 


2 jn(2n  + 1)  “i(ka) 
n=0  ^ 


[ Pn(cos  0 ’)Pn(cos  6 ) + 2 J -|P  “ p“  (cos  0 ')p“  (cos  0 ) cos  zn0'  ] 

n*-:  1 


a sin  0 'd0  'd0f 
(5. 133) 


Since 


277 

/ cos  m0'd0'  = 0 
0 


for 


m ^ 0 
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it  oo 


Nz  - -2a  / B2{  tlP2(t  j(cos  e •) 


00 


[ Z jn(2n  + 1)  (ka)P  (cos  8 *)P  (cos  8 ) ] sin  8 *d0 
n=0  n n n 


(5. 134) 


From  the  orthogonality  relationship  of  the  Legendre  functions, 


/ Pm  (cos  8 ) Pn(cos  8 ) (-Sin  8)  d8  = 0 


m / n 


2n  + 1 


m = n 

(5. 135) 


Therefore, 


00  W 

Nz  = a ^B2£+l^  ^2je+l(ka)  P2f+l(c0S  ( ITT^  + 3) 


4a  Jl  B2f  +1^  £ f2i+l(ka)  P2je+l(cos  9 ) 


(5.  136) 


Radiation  vector  N in  a spherical  coordinate  system 
is  related  to  the  Cartesian  components  by 


Nq  = (N  cos  0 + N sin  0)  cos  6 - N sin  6 
y x y “ 


N = -N  sin  0 + N cos  0 

0 xyy 


(5.  137) 


Therefore,  the  transverse  component  of  N becomes 
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N0  = TLcos  9 - N sin  6 
6 z 

N,  = 0 
0 


(5.  138) 


Due  to  the  symmetry  of  the  problem  with  respect  to  z 
axis,  it  was  anticipated  that  there  be  no  0 component  of  E fields,  at 
far-zone.  The  result  shown  in  equation  5.  137  and  5. 138  confirms 
the  physical  phenomena.  Since  the  direction  of  the  E fields  coincided 
with  that  of  the  N,  and  it  was  shown  that  = 0,  only  the  6 com- 
ponent of  E exists. 

Finally, 


(5. 139) 

[7|,cos  6 - Nz  sin  9]  ^ 


E = 
s 


jk7?Qe 


-jkR 


4jtR 


jMoe 


■jkR 


4ttR 


Equation  5.  139  is  the  expression  for  the  far-zone  electric  field  E 

s 

due  to  the  spherical  surface  current  that  was;  in  turn,  induced 
by  the  monopole  and  its  image. 
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5.  5. 1 Evaluation  ofTt.  From  Eq.  5. 131,  % is  defined  as 


H - / Z Z 


.n  2n+l 


o £=0 


2£+i  ^2£+rcosy  > 


jn(ka> 


P*(cos0’)  P*(cos0)  cos  O'  dO' 


since 


£m(x)  = (-l)m(l-x)m/2  — — l 
n dxm 


* d P (cos  O') 

Pn(cos0')  dS'  = d(cos9'), 


if  we  let  x = cost/',  the  integral  can  be  rewritten  as 


"ft-  ' 2a  •[,  £^0  B2t  + 1 P2£+l(x)  i"  ifsiy  i’n(ka) 


1 d Pn(x) 

P (cos0)  — 3 x dx 

n dx 


The  evaluation  of  this  integral  depends  upon  the  evaluation 


of  the  integral  I 


1 dPn(x) 

1 = /j  x P2£+l(x>  3x  * 


From  Ref.  17 


1 
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u 

d Pn(x) 

dx  - f2"'1)  pn.iW  + (2n-5)  Pn  3(x)  + (2n-9)  Pn_ g(x)  + . 

I 

(5. 142) 

1 

i 

also 

i 

(2n-l)xPnl(x)  = nPn(x)  + (n-l)Pu2(x) 

(5. 143) 

| 

From  (5. 142)  and  (5.  143) 

jS 

dPn(x) 

x dx  " fn  Pn(x)  + ^H'1)  Pn-2(x)^  + f(n“2)  pn_2(x)  + ("-3)  p 

n-4<x>l 

t 

i 

+ [("-■*)  pn.4(x)  + (n-5)  P„.5(x)?  + 

r> 

t 

- nPn(x)  + (2n-3)  Pn_2(x)  + (2n-7)  Pn_4(x)  + (2n-ll)  Pn  6(x)  + . . 

c 

<n/2 

= nP„<x)  + l (2n  - 4q  + 1)  P W 

q=l  n-zq 

(5.  144) 

» 

Thus 

? 

s> 

1 

T = [x  P2f +l(x>  [n  P„W  + (2"-3)  Pn  2(x)  + (2n- 7)  P^*)  + . . 

...  Jdx 

f 

(5.  145) 

5 

! 

From  the  orthogonal  relationship  of  the  legendre  polynomials, 

> 

l 

1 

pm(x)  Pn^x)  ^ = 0 for  n 

l 

i 

! 

. 2 

~ '2n+I  for  111  = « ■ 

(5.  135) 
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Using  the  orthogonality  shown  above,  the  integral  I can  be 


written  as 


1 1 <»/2 

I =/  npa+1(x)  Pn(x)  dx  + / Yj  (2n-4q+l)P2£+1(x)Pn_  2q«  * 


_ 26n-(2£+l)  /„  , n-21-1  2 

- (2{+1)  'm+ If+T  + \2n'  4 5 + 1 ) 2(22+1}  + 


2(2 '+1'  6n-<2£+J)  + 2 

qto  6n-2q-(2£+l) 


4£  +3 


(5. 146) 


where 


6 =1  for  m = 0 

m 

0 for  m/0  . 


Then, 


H-  2a  E B2Ulj  1 {21  + 1)(2£+2)  t2£  + l(ka^  P2Ul(c0Se) 


2(2f  + l) 
4C+3 


4a  B2J?  + 1 jn  ifn+T)  jn(ka)  Pi(cOS0) 


X 


..x2£  + l 1 


2a  f =0  Ba,+1<i)  * £ + 1 fat+l^  P22  + l(cos9* 


4a  {?0  ^ n?0  *”  i'0*’  P»(C0S(,)  • (5' 147) 
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Also  from  the  orthogonality  condition  of  the  second  integral  in  I ; i.  e. , 

n-  2*  - 1 
q = 2 


it  follows  that 


n = 2H  + 2q  + 1 


where 


i = 0,  1,  2,  ... 

q = 1,  2,  ... 

n = 2£+3,  2£+5,  . . . , cc  . 


Finally, 


oc 


- 2a  7 B i2£+1  p1  (cos  9) 

«L"  2 gQ  21+1 ] £ + 1 *2l+VC  S6) 


£^0  B2£+l  mJ+1  i2m+i  (2iil«ll)  t2m+l(ka)P2m+l(cos9) 


(5. 148) 


The  far- zone  electric  field  , due  to  the  surface  current  on  the  con- 
ducting sphere  induced  by  two  monopoles,  is: 
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6 


-jk  VQ  e 

snr 


-jkR 


ac°s9(£f  S^jj 


,24+1  fat+i(ka)  „i  , ,A 
m p2£+i<cos8) 


OC 


2 B 


1=0 


at+i 


V i2mtl  (4m+3)  4,  „ , „1  ,.\ 

m^+1J  (2m+l)(m+l)  d2m+l(ka)  P2m+l<cose>) 


2as,ne  £?0  B2£+i j2£+1  t2r.+l<ka>  P2£+l<cose>J 


which  becomes 


-jkR 


E9  = »R 
s 

B 


ka  cos  6 

2 — 


[•  B * 

Bi  ti(ka)  pj(cos  9)  " ~T  Bg(cos  9) 


+ ^ g(ka)  Pg(cos0) J 

[B1  (t  |3(ka)  P3(cos^  +TS  ^ka)  Bg(cos  0) 


ka  cos  t) 

+ 2 


15 

28 


• * 

^7(ka)  P7(cos0)  + | 


B3  (l5  t5(ka)  P5(cosf,)  - 15  j7M  P7(cos9)  + . . ^ 

B5  j7(ka)  p^(cose)  + ... 

• • 

B3  JVka^  ^(eos0)  + B5  ^g(ka)  P5(cos0)  - . . . 


^ + kasin#  jB,  ^(ka)  P£ 


(cos  9) 
(5.  149) 


The  ratio  of  the  first  term  to  second  term  is,  in  general, 
much  greater  than  1 as  was  shown  in  Tables  5.  1 and  5,  2 so  the  series 
converges  rapidly.  To  simplify  the  expression,  an  initial  term  is 
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taken  to  approximate  the  infinite  series.  It  will  be  shown  later  than  the 
subsequent  terms  do  not  contribute  significantly  toward  the  exact  value 
of  the  series. 

The  simplified  form  of  E is 


E. 


I n e 
m 'o 

»R 


-jkR 


(■ 


(ka)  Bj cos 0s in# 


) (3  " ^(5cos20-l)  ^(ka)^ 

(5. 150) 


where 


B' 


n 


2 B 
n 

^max 


The  numerical  calculation  of  the  far- zone  fields  E,  are 

6S 

performed  with  the  above  equation  (5. 150). 

5.  5.  2 Numerical  Evaluation  of  a Radiation  Pattern.  Since, 
in  Section  5. 3,  it  was  shown  that  the  far-zone  electric  field  due  to  two 
monopoles  separated  by  d = 2a  + h,  measured  between  center  to  center, 
has  an  expression  as  shown  in  equation  (5.  69);  i.  e. , 


E, 


i 7?  I e 
J 'o  m 


•jkR 


( 


cos(kd  cos  0)  - cos k(d-a)  cos(ka  cos  9)  + cos  0 sink(d-a)  sin(ka  cos  0) 

sin  0 


(5.69) 


The  total  electric  field  E at  far-zone  is 

V 


E*  = % * E»d 


(5. 151) 
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This  theoretical  evaluation  of  a far-zone  electric  field  due 
to  the  monopole  above  a semi- spherical  ground  plane  on  an  infinite 
ground  plane  was  conducted  mainly  to  evaluate  the  effect  of  the  detailed 
size  of  the  semi- spherical  ground  plane,  although  its  radius  was  much 
smaller  than  a wavelength.  Here,  a majority  of  experimental  studies 
were  performed  with  models  that  have  the  dimensions  of  a , corres- 
ponding to  the  radius  of  the  sphere,  less  than  or  equal  to  a quarter  of 
a wavelength  of  interest. 

The  expression  Efl  is  in  the  form  of  an  infinite  series 

where  the  coefficient  B and 

n 

metrical  size  of  the  ground  plane  and  a given  frequency.  Bn  is  a com- 
plex function  of  the  radius  of  the  sphere  and  the  antenna  height  in  terms 
of  wavelength;  i.e. , ka  and  kh. 


J^(ka)  are  constants  for  a given  geo- 


When  the  argument  of  a spherical  Bessel  function  is  much 
smaller  than  1,  it  can  be  evaluated  using  the  formula 


# 


n 

1 • 3 • 5 . . . (2n  + 1)  f0r  n = °’  2‘ • • 

z<<1  (5.152) 


Although  the  above  formula  is  for  z<<1,  the  ratio  of  these  functions 
between  n = 1,  and  n = 3 or  n = 5 are  fairly  close  to  the  actual  ratio 
of  these  functions  evaluated  with  exact  values  of  ^(z)  . 

Since  the  far-zone  electric  field  Eq  due  to  the  spherical 


surface  current  is  in  the  form  of  an  infinite  series  with  only  odd  n's  , 
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the  ratio  of  the  successive  terms  and  the  mature  of  its  convergence  is  impor- 
tant to  justify  taking  only  the  first  few  terms  in  the  actual  evaluation 
of  the  field  and  the  radiation  resistance. 

to  ^g(x)  are 

evaluated  by  using  the  approximate  formula  (5. 152),  and  by  using  the 
exact  values  for  x = 0. 5,  1.  and  1.  5.  Results  are  tabulated  in  Table  5. 1. 

Because  the  ground  plane  sizes  used  in  the  study  were 
limited  to  ka  < ^ which  corresponds  to  the  diameter  of  a semi- 
spherical  ground  plane  equal  to  or  less  than  one  half  of  a wavelength, 
the  ratios  only  of  those  spherical  Bessel  functions  with  an  argument  less 
than  or  equal  to  ^ are  of  interest  to  us. 


Ratios  of  the 


and^3(x) 


'"v^Arguments 
Ratio  N\ 

s. 

x = 0.  5 

x = 1.0 

x = 1.5 

/in 

Approx. 

144 

35.0 

15.55 

Exact 

138 

33.4 

14.0 

Approx. 

396 

99.0 

44.0 

' is(x>' 

Exact 

395 

97.2 

42.4 

Table  5. 1.  Ratios  of  ^(x)  . 
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Further  study  is  necessary*  to  determine  the  magnitude  of 
each  term  of  the  series.  Each  term,  in  addition  to  ^n(x) , contains 
Bn--  the  constants  evaluated  for  each  given  size  of  the  ground  plane 
radius  and  the  antenna  height.  Ratios  of  B^’  to  B^  and  B^  to  Bj!  can 
be  obtained  from  Table  2 for  several  typical  values  of  ka  and  kh. 


ka  = jt/2,  kh  = n/2 

ka  = t/3,  kh  = r/3 

ka  = x/6,  kh  = »/6 

V 

1.  298- j 0.742 

0.808  - j 0.078 

0. 342  - j 1. 30  x 10~3 

V 

-0.  247  + j 0.  738 

-0.057  +j  0.078 

0. 052  + j 1. 30  x 10~3 

V 

-0.046  + ] 4.  28  x 10"3 

0. 0287  + j 9.  66x  10-5 

0.031+j  1.  1 1 x 10"7 

Table  5.  2.  Coefficients  B ” 

n 

Finally,  the  behavior  of  Pn'cos  0)  and  P^^cos  6)  are 
graphically  shown  in  Figs.  5.  8 and  5. 9 for  0 <T  0 < 90°.  From 
Fig.  5.8,  the  magnitude  of  Pn(x)  is  to  be  always  less  than  or  equal 
to  1,  and  the  ratios  of  P^cos  9)  for  n = 1 and  3 or  n = 3 and  5 
are  always  less  than  6. 

11  can  be  concluded,  with  the  considerations  given  above, 
that  the  ratio  of  the  successive  terms  of  the  series  between  n = 1 
and  n = 3 is  always  less  than  0.  1 . Therefore,  it  is  jusitifed  using 
only  the  initial  term  from  each  infinite  series  in  Eq.  (5.  149)  or  (5.  153) 


P (cos  0) 
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0 (Degrees) 


Fig.  5.8.  Pn(cos  0)  versus  0. 
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V t,  v i2m+I  (4m  + 3)  1 v Di 
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(5. 153) 


00 


2 B2£+l  j 


i 21+1  : 


£=0 


^2£+l(ka)  P2£+1(C0S  G) 


Thus, 


j 77  I e 

c « 'o  m 

2iFR 


-jkR 


cos(kdcos6)  - cosk(d-a)  cos(kacos0)  + cosgsink(d-a)  sin(kacos0) 

sin  9 


] 


- J - 


f(ka)  BjMcos0  sine 


t 


_ 0 0 

j [^3  " \ <5cos20-  J)  ^ka)jj 


1(5. 154) 


The  plot  of  the  magnitude  of  the  term  inside  the  bracket  is 
a function  of  6 and  is  a radiation  pattern  for  a particular  a and  h. 

A few  of  the  patterns  are  calculated  and  are  plotted  below  in 
Fig.  5. 10  in  order  to  compare  with  the  experimental  results. 

Figure  5.  10  displays  a set  of  theoretical  radiation  patterns 
Eg  for  a monopole  and  for  its  image  in  the  upper  hemisphere.  This  is 
the  plot  of  Eq.  5.  154.  Since  the  radiation  patterns  are  symmetrical 
along  the  9 = 0°  axis,  only  the  right  half  of  the  complete  radiation 
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pattern  on  the  x - z plane  is  plotted  in  the  polar  coordinates. 

A minor  lobe  off  6=0°  axis  increases  with  the  size  of 
ground  diameter.  The  changes  in  the  ratio  of  the  intensity  of  the  major 
to  minor  lobe  is  shown  in  this  figure. 

There  are  a total  of  9 patterns  plotted  in  Fig.  5. 10.  Each 
row  is  a plot  with  a constant  distance  of  separation,  and  each  column 
with  a constant  frequency.  The  monopole  antenna  was  fixed  at  a quarter 
wavelength  for  the  highest  frequency  in  the  band.  Each  pattern  contains 
information  on  the  antenna  length  and  the  separation  distance  between 
the  two  antennas. 

Comparison  between  the  patterns  theoretically  obtained 
for  a monopole  and  its  image  separated  by  a diameter  of  the  spherical 
ground  plane,  but  without  the  presence  of  the  ground  plane,  to  those 
with  the  ground  plane  shows  that  adding  the  ground  plane  moved  the 
null  position  between  the  major  and  minor  lobes  toward  6 = 90°  axis 
at  higher  frequencies.  Also  at  lower  frequencies,  the  beam  width  was 
made  to  become  broader  by  adding  the  ground  plane.  Comparing  the 
theoretical  results  to  the  experime  al  results  also  show  that  the 
general  behaviours  of  radiation  patterns  are  similar  even  though  actual 
amplitude  do  not  seem  to  match  exactly,  'his  can  be  explained  by 
pointing  out  the  fact  that  the  initial  assumptions  made  concerning  the 
current  distributions  on  the  antenna  do  not  quite  agree  with  the  actual 
current  distributions.  Also  a finite  ground  plane  used  to  simulate  an 
infinite  ground  in  the  scale  models  caused  the  whole  pattern  to  be  shifted 
toward  0 = 0°  axis.  These  facts  are  theoretically  shown  by  Leitner  and 
Spence  (Ref.  3)  in  their  study  on  a monopole  with  a finite  ground  plane. 


^jt&rs<iZ>*' *r*5c 
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5. 6 Radiation  Resistance 

In  this  section,  the  radiation  resistance  is  calculated  as 
a function  of  the  antenna  length  and  grour.d  plane  diameter.  The  radi- 
ated power  is  calculated  using  the  far- zone  field  expression  that  was 
obtained  through  an  approximation  of  the  infinite  series.  The  far- 
zone  field  E0  was  given  in  Eq.  5.  154. 

The  convergence  of  Eq.5.  153  for  ka  < v , justifies 

using  only  the  first  few  terms  to  obtain  E^  . The  far-zone  electric 

u s 

field  due  to  a monopole  of  length  h and  a hemispherical  ground  plane 
of  radius  a is 


jtj  e 
'o 


-jkR 


2ttR 


ml 


cos(kd  cos  9)  - cos  k(d  - a)  cos(ka  cos  8) 
sin  0 


+ cos  0 sin  k(d  - a)  sin(ka  cos 
sin  0 ~~ 


(ka)  B j cos  0 sin  0 
2 


(ka)  - ^ (5  cos2  0 - 


(5.  154) 


The  first  and  second  terms  of  the  equation  in  the  brackets 
are  due  to  contributions  from  the  currents  on  the  two  monopoles  and  a 
spherical  ground  plane,  respectively. 
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Let 


F (q)  _ cos(kd  cos  0)  - cos  k(d  - a)  cos(ka  cos  9) 
1'  sin  6 


+ cos  Q sin  k(d  - a)  sin(ka 
sin  6 


and 


F2(6}  = 


(ka)  cos  9 sin  9 [ ^(ka)  +^(5  cos*  9 - 1)  ^3(ka)j 


Then,  the  Eq.  5. 154  can  be  written  as 


% *'jkR  r 

E8  ■ —m  lm  [Fl<«>  + J B 


1 F2(8)] 


Letting  = a ^ + j/3  , Eg  can  be  written  as 


-jkR 

Jr?oe 

E0  " " 2irR  lm 


[ Fj(9)  + (a  j + j(3j)  F2(9)l 


-jkR 


2ttR 


» 

'm  [ [F,<»> +/»,*!<«]*  J“l  F2(9 


The  far- zone  magnetic  field  intensity  is,  from  H = — (SxE) 

n / 

'o 


H. 


_ J e 


■jkR 


277  R 


m 


[[Fl(#)+ W»]-  J«j  *,<•>] 


cos  9) 

(5.  71) 


(5.  155) 


(5.  156) 


(5.  157) 


(5.  158) 
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jkR 


2trR 


m 


[[Fl<8>  + h F2<e>] 


+ j al  F2(6) 


j 


Radiation  resistance  R is  defined  as 

r 


(5.  159) 


R = — L_  (g  r i (E  x h*)  • dS 


I 


2 y «e  2 


rms 


~ i #ReKv|R8  sln0ded 

rms  * ' 


1 - M 


2 I 2 w 4ir2  R2 
rms  ™ n 


R2sin  9 dO  d0 


^ f [ Fl2  <9>  + 2<31  F!<»)  Fj,(9)  + (3j2  F/  (9)  + 0l2  F/  (9)1 


sin  0 d9  d0 


r?0  tt 

= 2iT  / [F12(9)+  2/Sj  F1(9)  F2(9)  + (Oj2  + ) F^  (9)J  sin  9 d9 

(5.  160) 

Since  ri0  = 1120  v in  free-space, 

Rr  = 60  / [fi2(0)J  sin  9 d0  + 120^/ F^O)  F2(9)  sin  9 d0 

+ 60(ax2  + ) j F22  (0)  sin  0 d9 

(5.  161) 


s*b 
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Let 


where 


R — R + R + R 
r rl  r2  r3 


ir 


Rr  =60  / F^G)  sinO  d0 


7T 


Rr  = 120  /3X  / Fj(9)  F2(0)  sin  9 dG 
2 o 


Rr  = 60(Ol2  +0t8)  f F22  (G)  sin  G dG 
3 o 


5.6.1  Evaluation  of  R . From  Eq.  FJ9)  is  given  as 

1 1 


F ,q\  cos(kd  cos  Q)  - cos  k(d  - a)  :os(ka  cos  9) 
1'  sin  0 


+ cos  9 sin  k(d  - a)  sin(ka  cos  0) 
sin  9 


(5.  162) 


(5.  163) 


(5.71) 


Therefore, 
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Let  ju  = cos  9 , d/i  = -sin  6 d0  and  d = a + h and 

solve  integrals  from  (l)  through  © separately.  Then  sum  of  the  results 

will  give  R . In  evaluating  each  integral,  logarithmic  singularities 
rl 

appeal  These  singularities  will  later  cancel.  However,  in  order  to 
show  the  convergence  of  the  integrals,  the  limits  of  the  integration  have 
been  changed  from  /u  = - 1 and  fi  - 1 to  (i  = -1+e  and  /i  = 1 - e 
and  then  € was  made  vanishingly  small. 

Then, 

(T)  . lim  /£  C°s2  Ma-^^  d, 

€ - O - l+  c 1 ' M 

Since 

TH?  2 + m 1 - **  ) 

and 

cos2  k(a  + h)  p = Y 

©=  i [rnr  t^t]  ( 1 + cos  2k(a  + h)  ")dfI 

When 

F(jm)  = F(-m) 


( 


1 + cos  2k(a  + h)  p, 


) 
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Therefore, 


/ fcm)_  . Jm 

a 1 + 11  i1-" 


© 


D = 2 lim  f 

€ — O - l+£ 


1 + cos  2k(a  + h)/jL 

1 + M 


d m 


1 ,imr  / <!*_  + f cos  2k(a  + h)  M 'd1 

2 f-oL-l+(  1 + (I  -1«  I + " J 

After  the  change  of  variables,  1 + /i  = u for  the  first  integral 
and  1 + n = u and  2k(a  + h)u  = v for  the  second  integral.  (T)  can  be 
written  as 


1 

“ 2 

lim 

/ 

du 

— + cos 
u 

2k(a  + 

e -o 

L e 

+ sit 

\ 2k(a 

+ h) 

4kja+h) 

sin  v 

if 

2k(a+h)e 

V 

Ci 

(4k  (a 

+ h)Y 

• lim  Ci 

(2k(a  + 

\ 

/ 

c-  0 

\ 

dv 


dv 


1 

2 


f n 2 - lim  ( n f + cos  2k(a  + h) 
e-  o 


+ sin  2k(a  + h)  Si  [4k(a  + h) 


i ^4k(a  + h)^  | 
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I 

;■ 

i 

( 

n 


I 


i 


where 

x 

Si(x)  = f -1-1— - dv  = sine  integral 

6 v 

x 

f COS  V 

Ci(x)  = dv  = cosine  integral 

*'  v 
x 

In  the  solution  for  (T)  , singularities  associated  with  logarithm 
and  cosine  integral  are  left  as  they  are  and  it  will  be  shown  later  that  they 
cancel. 


Using  a similar  technique,  it  can  be  shown  that 


( 2 ) = cos2  (kh)  lim  f dp 

^ wb  i-m4* 


^ cos2  (kh)  | t n 2 - lim  fn  c + cos  2ka  Jci(4ka)  - lim  Ci(2ka  e)j 
( c —o  e—o 


+ sin  2 ka  Si(4  ka) 


) 


1 


7i*VJHR 


X~€ 

(sin*  kh)  lim  f 

€— *0  t 
-l+€ 


M2  sin*  (ka  m)  da 
I - m2 


adding  and  subtracting  sin2(kaju)  in  the  numerator. 
= (sin*  kh)  lim  J -*n*  (ka  jj)  “ U - M2 ) sin*  ka  u 


C-0  -1t€ 


1-  M‘ 


= (sin2  kh)  lim  [ f -g‘-?a  **£  dll  -Y  sin2  ka  n dJ 

<-oL-l«  1-M  -i«  J 


The  result  of  the  integration  is 


® = 1 (sin*  kh)  £n  2 - lim  i n c - 2 (l 

€-0  V 


sin  2 ka ' 
2 ka  j 


- cos  2 ka  JCi(4  ka)  - lim  Ci(2  ka  e ) - sin  2 ka  Si (4 ka) J 

L €-0  J 


® = -2  cos(kh)  lim  jf  g_os  k(a  + h)  cos(ka  u) 


c—*o  -l+e 


using  an  identity  that 


1-  M* 


cos  A cos  B = ^ j^cos(A  + B)  + cos  (A  - B)J 


180 


(T)  = - cos(kh)  lim 


e— *o 


f cos  kj_2a  ±h)jx  d f cos  khji_ 
-1«  1 + 11  -3«  1 + ,J 


cos(kh) 


cos  k(2a  + h)  |^Ci^2k(2a  + h^-  lim  Ci  ^k(2a  + h)  e jj 

] 


^2k(2a  + h)J  + cos  kh  Ci 


+ sin  k(2a  + h)  Si  (2k(2a  + h)j  + cos  kh  Ci(2  kh)  - lim  Ci(kh  e) 


+ sin  kh  Si  (2  kh) 


(5)  = 2 sin  kh  lim  f 
e— o - 1+c 


l±  cos  k(a  + h)  i±  sin  ka  m d/i 
1-  M2 


using  the  identities 


1 - M 


2 2 


L_) 

\1  - M 1 + M / 


and  coc  A sin  B = ||sin(A  + B)  - sin(A  - B)] 


and  also  f dfj.  = - f d/i  when  F(ju)  = -F(-| 1) 

a ^ a + ^ 


the  result  of  the  integral  (?)  becomes 


(F)  = sin  kh 


sin  k(2a  + h) 


jci^2k(2a  + h)j  - lim  Ci  (k{2a.  + h) 


- cos  k(2a  + h)  Si  (2k(2a  + h) 


^2k(2a  + h)^  - sin  kh  Ci 


Ci(2kh)  - lim  Ci(kh  e) 
e-*o 


+ cos  kh  Si(2kh) 
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Similarly, 


1-e 


® = -2  cos  kh  sin  kh  lim  / cos  te  H sin  to  M 


d/i 


€-*0  - 1+€ 


1-  M‘ 


g-  sin  2 kh 


fc 


cos  2 ka  Si(4  ka)  - sin  2 ka  Ci(4  ka)  - lim  Ci(kh  e ) 

L €-0 


1 


Cosine  integral  Ci(x)  can  be  expanded  into  a power  series 
such  as  (Ref.  18,  p.  232). 


oo  ( ,>n  „2n 

Ci(z)  = y + f n Z + £ _2n(2n)! 


(5.  165) 


If  argument  of  the  cosine  integral,  z,  is  small 


lim  Ci(z)  ~ y+  fnz 
z — o 


(5.  166) 


where 


y = 0. 577  = Euler's  constant 


Therefore,  lim  Ci(2k  ae),  for  instance,  can  be  broken  up  into  two  terms; 
e-  o 

i.  e. , 


lim  Ci(2  ka  e)  « y + £ n{2  ka)  + lim  fn® 
e -o  e — o 


(5.  167) 
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Applying  the  technique  shown  inEq.5.166and  adding  the  results (T) 
through  (tT),  it  can  be  shown  (Appendix  E)  that 


30 


2 y + 2£n(2  kh)  - 2 Ci(2  kh)  - 2 sin2  kh  (1  - 


sin  2 ka 
2 ka 


+ cos  2k(a  + h)  |^Ci(4  ka)  + Ci^4k(a  + h)^  +2  £n  k(2a  + h)  - £n(2  ka) 

i |^2k(2a  + h)  J j + sin  2k(a+  h)  ^Si(4  ka)  + Si^k  (a  + h)^ 


- £n  2k(a  + h)  - 2 Ci 

- 2 Si^2k(a  + h)j 


(5. 168) 


If  we  impose  a condition  that  a = 0 and  h = m ^ where 
m is  an  odd  integer,  the  result  matches  with  that  of  Stratton's  (Ref.  14, 
p.  444)  where  he  found  a radiation  resistance  of  a linear  dipole  at  reson- 
ance. 


When 


a 


0 and  h = ^ 


Eq.  5.  168  reduces  to 


30 


£ y + £ n 2n  - 


(5.  169) 


which  is  the  exact  replica  of  Stratton's  expression 

R = 30  J^£n  277  y*  - Ci  2m  irj  with  m = 1 (5.  170) 

y*  = 1.7811 

It  can  further  be  shown  that  when  h = 0,  R =0.  This, 

rl 

of  course,  must  be  true  for  physical  reasons  and  this  gives  another 
confirmation  of  the  correctness  of  this  theoretical  result. 


183 


5.6 . 2 Evaluation  of  R 


From  Eq.  5. 163, 


Rr  = 120^  / F1(0)F2(0)  sine  d9 
2 o 


_ ^20  3 f cos(^d  c°s  6)  - cos  k(d  - a)  cos(ka  cos  Q) 

P1  " [ sin  6 


+ cos  £ sin  k(d  - a)  sin(ka  cos  Q)  1 


sin  6 


Tka  cos  0 sin  0 {3  ^(ka)  - 1 (5  cos2  Q - 1)  ^3(ka)jl 

*■  2 J su 


G(ka)--y-[^  (ka)  - -3^  (ka)  j 

H(ka)  = (y)  (y  (ka))  = (ka)  j3 (ka) 


FgW  = G(ka)  sin  6 cos  0 + H(ka)  sin  9 cos3  0 


R can  now  be  expressed  as 
2 


% “ i2°  \ / |^cos  kd  x " cos  k(d  - a)  cos  ka  x + x sin  k(d  - a)  sin(ka  x)J 
[ G(ka)  x + H(ka)  x3]  dx 

Every  term  in  the  above  integral  is  an  odd  function  of  x with 
respect  to  x = 0 . Therefore,  the  result  of  integration  becomes  zero. 

Therefore 


R s 0 
r2 


(5. 174) 


5.  6.  3 Evaluation  of  R . 

r3 


{r  =60(0^+^)  / F2(O)sin0d© 

3 n 


60(al2  + /Sj2)  / J^G(ka)  sin  0 cos  0 + H(ka)  sin  0 cos3e]‘ 


sin  0 d0 


v r 

= / [G2(ka)sin20cos20  + 2G(ka)H(ka)sin2Ocos40 


H2  (ka)  sin*0  cos60j  sin  0 d9 


(5.  17?) 
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Letting 

cos  0 = x -sin  9 dO  = dx 

sin2  6 = 1 - x2 

Rr  = 60(a12  + 0t2 ) f [g2  (ka)(l  - x2 ) x2  + 2G(ka)  H(ka)(l  - x2)  x4 
3 -1 

+ H2  (ka)(l  - x2  ) x6]  dx 

= 60(a12  + 0^ ) f G2  (ka)  x2  + [ 2G(ka)  H(ka)  - G2  (ka)]  x4 
+ f H2  (ka)  - 2G(ka)  H(ka)]  x6  - H2  (ka)  x8  dx 
= 60 (a  j2  + 0 2 ) G2  (ka)  (|-)  +[2G(ka)  H(ka)  - G2  (ka)]  ( 

+ [h2  (ka)  - 2G(ka)  H(ka)](|)  - H2  (ka)  (|) 

= 60 (a  2 + 0 x2 ) — - G2  (ka)  + ^ G(ka)  H(ka)  + ^ H2  (ka) 

G2  (ka)  2G(ka)  H(ka)  H2  (ka) 

15  + 35  'r  63 


= 240(ttl2  +/312) 


(5.  176) 


cn|  to 


Finally,  ihe  total  radiation  resistance  of  the  system  with 
monopole  and  a hemispherical  ground  plane  is  obtained  by  summing 

Rrj  9 Rr2  and  Rr3  * However»  this  sum  represents  a radiation 
resistance  of  the  actual  antenna  system  and  its  image.  Because  only 

a half  of  the  power  evaluated  previously  is  actually  radiated,  mono- 
pole and  hemispherical  ground  plane  radiating  into  a half  space,  the 
final  radiation  resistance  is 


Rad.  Resistance  = 4(R  + R +R  ) 

AT  r'  ' 


1 r2  r3 


= 15  2 y + 2 in(2  kh)  - 2 Ci(2kh)  - 2 sin2  kh  (1  - -gin  2 ka) 


+ cos  2k<a  + h>  jci(4  ka)  + Ci^4k(a  + h))  + 2 in  k(2a  + h) 

- I n(2  ka)  - £n  2k(a  + h)  - 2 Ci  [2k(2a  + h)]J 

+ sin  2k(a  + h)  j^Si(4  ka)  + Si^4k(a  + h)J  - 2 Si  ^2k(a  7h)j}J  j 

+•  120(0*  2 +0  *)  ( , 2G(ka)  H(ka)  H2(ka)| 

1 1 / 15  35  + 63  1 


+ sin 


From  Eq.  5.  177,  radiation  resistances  of  several  different 


combinations  of  ka  and  kh 


are  numerically  calculated  and  tabulated. 
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The  result  of  numerical  evaluation  of  Eq.  5. 177  is  given  in 
a ,raph?cal  form  in  Fig.  5. 11,  where  radiation  resistance  is  plotted  as 
a function  of  frequency.  Each  curve  represents  a different  size  of  semi- 
spherical  ground  plane.  When  the  radius  of  a ground  plane  is  zero, 
the  resulting  radiation  resistance  corresponds  to  that  of  a monopole  on 
an  infinite  ground  plane. 

Results  shown  in  Fig.  5. 11  indicates  that  there  exists 
definitely  a peaking  effect  on  the  radiation  resistance  as  the  radius 
of  the  ground  plane  is  changed.  Comparing  with  the  monopole  resonating 
at  30  MHz  over  an  infinite  ground  plane  the  radiation  resistance  becomes 
larger  with  a semi- spherical  ground  plane  well  below  the  resonant 
frequency.  The  peaking  seems  to  occur  at  ka  = 1 and  this  conclusion 
has  been  drawn  mainly  from  the  results  of  numerical  calculations. 

Fig.  5.  12  and  5.  13  compare  these  theoretical  results  with 
experimentally  obtained  input  resistances.  Of  course,  we  are  not 
comparing  the  same  resistances,  namely  the  radiation  resistances. 
However,  assuming  that  the  loss  is  small,  the  input  resistance  should 
be  similar  to  the  radiation  resistances.  Some  of  the  discrepancies 
shown  in  this  comparison  can  also  be  explained  with  the  discrepancies 
in  the  assumed  current  and  the  actual  current  on  the  antenna.  However, 
the  peaking  effect  is  shown  to  exist  using  a small  semi- spherical  ground 
plane  with  a monopole. 
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Frequency  (MHz) 


Fig.  5.  12.  Theoretical  radiation  resistance  and  experimental 

input  resistance  for  a monopole  with  a hemispherical 
ground  plane 
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Frequency  (MHz) 


Fig.  5. 13.  Theoretical  radiation  resistance  and  experimental 

input  resistance  for  a monopole  with  a hemispherical 
ground  plane 
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Since  radiation  resistance  is  defined  as  the  ratio  of  radiated 
power  to  a square  of  the  input  current  at  the  terminal  where  it  is  mea- 
sured, the  current  must  be  that  current  at  the  terminal  and  resulting 
radiation  resistance  will  be  the  radiation  resistance  at  that  point. 

At  the  beginning  of  the  theoretical  analysis  for  this  problem, 
it  w as  assumed  that  the  current  distribution  on  the  antenna  itself  is  in 
the  sinusoidal  form. 

From  Eq.  5.  59  , the  current  on  the  antenna  is  given  as 

I(z)  = I sin  k(d  - z) 
max 

where  at  z - a , .An  k(d  - z)  = 1 at  30  MHz. 

Therefore.,  the  input  currents  at  the  base  of  the  antenna 

(z  = a)  would  be  different  from  I for  frequencies  lower  than  30  MHz. 

max 

We  must  use  these  currents  to  evaluate  the  radiation  resistance  at  the 
input  terminal  for  each  iroqvenoy. 

It  seems  as  though  one  nan  have  a radiation  resistance  as 
large  as  desired  by  making  the  feedpoint  located  at  a point  where  current 
is  approximately  zero.  One  example  would  be  a full  wave  dipole  which 
has  a zero  input  current  at  the  center  of  the  dipole.  Theoretically, 
then,  it  has  an  infinitely  large  rac  iation  resistance.  However,  in 
practice  an  actual  antenna  is  not  infinitesimally  thin,  which  was  the 
assumption  used  for  deriving  a theoretical  result,  and  the  current  at 
a minimum  point  is  not  zero.  Nevertheless,  the  radiation  resistance 


| at  a current  minimum  may  in  practice  be  very  large,  possibly 

I thousands  of  ohms. 

f 

In  view  of  the  above  argument,  the  radiation  resistance 
cannot  be  used  as  a sole  measure  of  an  antenna  performance.  The 
reactive  component  of  an  input  impedance  is  another  factor  which  should 
be  considered.  When  an  input  current  is  fixed  at  a certain  value,  both 
the  resistive  and  reactive  components  of  the  input  impedance  determine 
the  voltage  required  to  maintain  the  current.  Therefore,  it  is  necessary 
to  consider  a ratio  of  reactive  component  and  resistive  components, 
or  the  Q of  the  antenna  in  order  to  decide  whether  increases  in  radiation 
resistance  actually  increase  the  antenna  efficiency.  Input  current  and 
voltage  determines  the  input  power  and  radiation  resistance  determines 
the  radiated  power.  The  ratio  of  these  quantities  determines  the  efficiency 
of  the  antenna. 

In  addition  to  the  efficiency  considerations,  the  Q of  the 
antenna  system  determines  broadband  tuning  possibilities  of  the  antenna. 
When  the  antenna  is  electrically  short,  in  particular,  the  reactive 
component  is  usually  very  large  and  the  resistive  component  is  very 
small  giving  very  high  Q values.  When  the  reactive  component  is  large, 
a wideband  tuning  becomes  very  difficult,  because  one  has  to  find  a 
way  for  a wide  range  variation  in  the  conjugate  reactive  component  to 
cancel  out  the  large  radiative  reactances.  Also,  a small  value  of  resistive 


component  will  cause  a difficulty  in  transforming  up  to  a characteristic 
impedance  of  transmission  lines  normally  available. 
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The  ratio  of  a capacitive  reactance  to  a resistance  of  the 
antenna  is  defined  as  Q,  for  a monopoie  above  an  infinite  ground 
plane,  Q of  the  antenna  increases  rapidly  as  frequencies  go  down 
below  resonance.  Typically,  for  a monopole  one  quarter  wavelength 
long  at  30  MHz,  with  the  antenna  or  diameter  ratio  of  approximately 
400,  Q at  5MHz  becomes  1400  and  gradually  reduces  to  0 at  30MHz. 
The  resistive  component  increases  from  2 ohms  at  5 MHz  to  36  ohms 
at  30Mhz.  However,  from  the  experimentally  obtained  data,  Q of 
a monopole  one  quarter  wavelength  long  at  30  MHz  placed  above 
a finite  circular  disc  ground  plane  of  2.  5 meters  in  diameter  which 
is  placed  1.  25  meters  above  an  infinite  natural  ground  is  shown  to 
vary  between  43  to  0.  It  is  mainly  due  to  smaller  reactive  components 
and  large  resistive  components  over  the  frequency  range  of  SMHz 
through  30MHz  for  a monopole  over  a finite  ground  plane. 

The  resistive  component  of  the  input  impedance  contains 
loss  due  to  several  sources  in  addition  to  radiation  resistance.  In 
order  to  make  an  educated  judgment  on  the  variation  of  the  resistive 
component  as  a function  of  frequency  whether  it  can  be  attributed  to 
an  increase  in  radiation  resistance  or  not,  a theoretical  result  of 
a radiation  resistance  for  a mon  jole  over  semi-  spherical  ground 
planes  has  been  obtained.  Using  these  theoretical  results  of  radiation, 
resistance  and  the  experimentally  obtained  reactive  components,  Q 
of  the  antenna  has  been  obtained  between  5MHz  and  30MHz.  It  shows 
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that  the  Q of  the  antenna  varies  between  25  to  1 with  X at 

V 

10 MHz  of  850  ohms  compared  with  900  ohms  for  a monopole  above 
an  infinite  ground  plane. 

A conclusion  can  be  drawn  from  the  above  studies  that 
by  using  a monopole  with  a finite  ground  plane  above  an  infinitely 
large  ground  plane  , a significant  advantage  can  be  realized  in  rad- 
iation resistance  as  well  as  a broad  band  matching  at  the  input  terminals. 
Through  this  study,  a qualitative  conclusion  of  an  antenna  behavior 
with  two  ground  planes  has  been  possible.  However,  an  exact 
behavior  of  a monopole  with  a disc  ground  plane  above  an  infinite 
ground  can  only  be  predicted  with  a new  mathematical  model. 


CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 
6.  1 Conclusion 

One  of  the  major  objectives  of  this  study  was  to  determine 
how  a finite  ground  plane  used  in  addition  to  an  infinite  lossly  ground 
affects  performa  nee  of  a monopole  antenna.  Attention  was  mainly 
directed  toward  the  effect  of  ground  planes  equal  to,  or  smaller  than, 
a half  wavelength  in  diameter.  If  there  is  to  be  any  advantage  in 
using  these  finite  ground  planes,  they  cannot  be  physically  bulky.  At 
frequencies  5MHz  to  30  MHz,  where  the  study  was  performed,  a half 
wavelength  at  the  upper  end  of  the  frequency  band  is  equal  to  5 meters. 

In  studying  this  antenna  system,  the  length  of  monopole 
was  fUed  at  2.  5 meters,  which  is  a quarter  wavelength  at  30  MHz. 
Impedance  measurements  show  that  the  real  part  is  a strong  function 
of  both  the  ground  plane  diameter  and  its  location.  Also  it  shows  a 
peaking  effect  at  particular  frequencies,  which  depends  upon  the 
diameter  of  the  ground  plane  and  its  location  above  the  infinite  ground. 

Using  a theoretical  model  with  hemispherical  ground 
planes  suostituted  for  flat  discs,  it  was  established  that  the  sharp 
increase  in  resistive  component  was  largely  due  to  an  increase  in  radia- 
tion resistance.  From  the  theory,  it  seems  that  at  least  one  peak 
exists  at  frequencies  where  ka  - 1.  It  was  not  possible  to  draw  a 
clear  conclusion  of  this  nature  for  an  antenna  with  a finite  disc  ground 
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plane.  When  ka  « 1,  the  distance  from  the  base  of  the  antenna  to  an 
infinite  ground  a eng  the  path  of  the  spherical  surface  is  a quarter 
wavelength.  This  provides  a voltage  maximum  or  current  minimum 
point  at  the  base  of  the  antenna,  regarding  the  infinite  ground  plane 
as  a zero  potential  surface.  Therefore,  a larger  input  impedance  or 
larger  resistive  component  is  realized  at  ka  = 1.  However,  based 
on  the  experimental  results,  a general  conclusion  that  a peaking  effect 
occurs  at  frequencies  approximately  determined  by  + a^  C where 
C is  a function  of  the  ratio  It  can  also  be  concluded  that  the  Q 

of  the  antenna  is  generally  much  lower  than  that  of  the  same  antenna 
over  an  infinite  ground  plane,  providing  some  advantage  in  designing 
a tuning  network  at  frequencies  well  below  resonance. 

Curves  showing  the  input  reactance  with  a finite  disc- 
ground  plane  of  various  sizes  exhibit  a region  where  the  reactance 
variation  as  a function  of  frequency  is  remarkably  small.  This 
phenomenon  is  usually  associated  with  a resistive  maximum.  Clearly, 
this  type  of  frequency  response  is  very  much  better,  from  the  point 
of  view  of  broad-band  operation,  than  the  response  of  either  a half- 
wave or  a full- wave  dipole.  With  an  appropriately  designed  reactive 
matching  network  the  standing  wave  ratio  on  a line  terminated  in 
the  antenna  with  a finite  ground  plane  over  an  infinite  ground  can  be 
minimized  over  a wide  frequency  range. 

Because  of  the  assumption  made  in  the  theoretical  studies 
that  the  current  distribution  on  the  monopole  is  in  sinusoidal  form 
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and  independent  of  the  surface  current  on  the  ground  plane,  actual 
current  measurements  were  performed  using  scale  models.  Current 
distribution  at  the  higher  end  of  the  frequency  band  approximated 
closely  the  initial  assumption  of  a sinusoidal  distribution.  However, 
as  frequency  is  decreased,  particularly  for  the  monopole  with  a hemi- 
spherical ground,  the  actual  current  deviates  appreciably  from  the 
original  assumption.  At  the  lower  end  of  the  frequency  band,  the 
current  measurement  shows  almost  a constant  amplitude  similar  to 
Hertzian  dipole. 

The  theoretical  results  obtained  here  for  the  radiation 
resistance  assumed  that  the  current  distribution  had  a triangular 
form  at  low  frequencies.  These  results  would  be  improved  by 
assumption,  in  the  theory,  of  a current  distribution  more  closely 
approximating  that  observed  experimentally.  This  effect  is  much 
more  noticeable  at  the  lower  end  of  frequency  band.  Experimentally 
measured  input  resistances  confirm  these  facts.  In  Chapter  5,  a 
comparison  is  made  between  the  theoretical  and  experimental  radia- 
tion and  input  resistances.  The  experimental  values  are  always 
larger  than  the  theoretical  values. 

From  the  results  of  a theoretical  study  using  a hemi- 
spherical ground  plane  and  the  experimental  results  obtained  for  a 
monopole  with  a finite  disc  ground  plane  above  an  infinite  natural 
ground,  it  can  be  concluded  that  a marked  increase  in  radiation 
resistance  results  with  a finite  disc  ground  plane,  below  the  resonant 
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frequency. 

The  major  effect  of  the  finite  ground  plane  upon  the  antenna 
radiation  pattern  is  an  emergence  of  side  lobes  when  the  ground  plane 
diameter  becomes  appreciable  compared  with  a wavelength.  In  general, 
the  appearance  of  a side  lobe  is  not  desirable,  because  the 
power  radiated  is  actually  deviated  from  the  main  lobe  where  it  should 
be  concentrated.  Also  the  beam  width  measured  between  3db  points 
from  the  position  of  the  maximum  amplitude  is  generally  broader 
than  that  of  an  antenna  on  an  infinite  ground  plane.  Theoretically,  it 
has  been  found  that  a spherical  surface  current  makes  only  a very 
slight  contribution  to  the  far- zone  electromagnetic  field. 

In  summarizing  the  results  of  this  study,  it  is  concluded 
that  definitely  improved  performance  can  be  obtained  with  an  electri- 
cally-short  monopole  antenna  by  operating  it  with  a physically  small 
finite  ground  plane  above  natural  ground.  This  improvement  occurs 
in  the  radiation  resistance  and  with  respect  to  matching  network 
design  for  a broad- band  operation.  These  results  have  bearing  on 
the  design  of  short  monopole  antennas  for  operation  of  ground- 
based  vehicles,  in  which  case  the  antenna  may  be  considered  to  be 
mounted  on  a small  ground  plane  above  natural  ground. 

G.  2 Recommendation  for  Future  Work 

First,  in  the  analysis  of  this  problem,  the  infinite  ground 
below  the  finite  hemispherical  ground  was  assumed  to  be  perfectly 
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conducting,  permitting  use  of  a conventional  image  technique  for 
solving  the  problem.  The  loss  due  to  the  finite  conductivity  of  the 
natural  ground  has  also  been  neglected.  However,  in  order  to  isolate 
an  exact  radiation  resistance  from  the  measured  impedance  shown 
in  Chapter  2,  it  is  necessary  to  take  the  loss  due  to  the  natural  ground 
into  account. 

Also,  a new  theoretical  model  consisting  of  an  antenna  on 
a small  ground  plane  above  a larger  ground  plane,  could  possibly  give 
more  accurate  results  than  have  been  obtained  here  through  assumption 
of  a hemispherical  ground  plane. 

Second,  since  the  current  distribution  measured  shows 
deviation  from  the  original  sinusoidal  form,  the  measured  current 
may  be  taken  as  a given  current  distribution  and  may  be  decomposed 
into  Fourier  Series  to  add  an  effect  of  each  term.  In  this  way,  the 
experimental  result  may  be  made  to  match  better  against  theoretical 
results. 

Third,  because  the  ultimate  application  of  this  analysis 
is  for  a vehicular  mounted  antenna  where  a body  of  a vehicle  can  be 
taken  as  an  equivalent  ground  plane  of  a finite  diameter,  it  is  necessary 
to  study  further  to  determine  how  a particular  size  of  a vehicular  body 
and  shape  and  location  of  the  antenna  correspond  with  an  antenna  on  a 
finite  disc  or  semispherical  ground  plane  above  an  infinite  natural 
gound. 

Finally,  a detailed  study  of  the  measured  input  impedance 
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given  in  Chapter  2 and  deriving  equivalent  circuits  for  each  case  may 
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APPENDIX  A 
RECIPROCITY  THEOREM 

Let  a source  or  a distribution  of  sources  maintain  a current 
density  with  associated  electromagnetic  fields  E^  and  . A 
second  source  maintains  a current  density  J 2 with  associated  fields 
Eg  and  Hg  . Then,  according  to  Lorentz's  reciprocal  theorem: 

V • (Ej  x H2)  - V • (Eg  x Hj)  = Jx  • Eg  - J2  • Ej  (A.  1) 

Let's  assume  that  the  two  sets  of  fields  are  defined  in  free 
space  and  are  bounded  by  the  same  geometrical  surfaces.  Applying 
volume  integration  on  both  sides  of  Eq.  A.  1 with  dV  as  volume  element, 
the  integration  performed  over  a volume  V in  empty  space  that  is 
bounded  by  a closed  surface  S can  be  written  as 

Iff  V ■ (f  J x H2)  dV  = ///  V • (Eg  x Hj)  dV  (A.  2) 


where 


///(Jj-  Eg-  J2-  Ej)dV  = 0 (A.  3) 

because  all  volume  densities  of  current  vanish  in  empty  space. 

From  the  divergence  theorem 
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///  V * (Ej  x H2)  - (E2  x Hj)  dV  = ff  x H2  - E2  x Hj)  * dS  = 0 

(A.  4) 

Where  dS  = n dS  and  n is  an  external  unit  normal  vector 
of  the  surface  toward  volume  V . The  only  restrictions  of  the  two  sets 
of  fields  (E^  , H^)  and  (Eg  , Hg)  are  that  they  satisfy  Maxwell’s  equations, 
the  volume  integration  is  performed  in  empty  space  V and  the  boundary 
surfaces  are  in  the  same  geometrical  shape  for  both  sets,  but  not  neces- 
sarily the  same  physical  properties. 


APPENDIX  B 


PROOF  AND  DERIVATION  OF  EQ.  5.  56 


Equation  20  was  written  as 


E(x.  y,  z)  = 


At  far-zone, 


-T-  J f [jo? /x(n  x H)  0 + ( n x E)  x V * 0 + (n  • E)  V ’ d]  dS* 
w g 

(5.  50) 


r 


/ 


(x  - x’)2  + (y  - y’)2  + (z  - z’); 


IR-R'I  (5.39) 


can  be  simplified  as 


r = R - 


R • R’ 
R 


= R - R • R’ 


Therefore, 


0 


gjkR  - jk$  • R’ 
_ 


(5.  55) 


when 


I R I »IR'I 
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Then 


V'0  * 


aikR 

"R 


. -jkeikR 


R 


e-jlcR  R'  A 


(5.47) 


where  the  prime  notation  on  operator  V indicates  that  V is  operated 
on  the  primed  coordinate  system. 

In  Eq.  5.  50,  the  radial  component  of  E (x,  y,  z)  can  be 

written  as 


Er(x,  y,  z)  = 


1 

~w 


//s 

s 


jwjLt^x  H)  - jk(^ 


E)FT 


JkR 


jkR  • R1 


dS’ 


(B.l) 


Consider  the  following  integral 


§ V'  x|He”ikR‘  R’  y ds'  = 0 (B.  2) 

Where  S is  a closed  surface  bounding  a volume  V . 

Since 


Vx(0a)  = V0xa  + 0Vxa 


V’  x 


77  -jk^ 

xHte 


R 


V'  x H 


(B.3) 
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Therefore,  Eq.  B.  2 becomes 


Jf  fn  • (V’  x H)e"ikR  ' R 


o 


(B.4) 

where 


dS  = n dS 


Knowing 


and 


V’e 


-jlcR  4 R* 


A 


R 


“n-  (RxH)  = -R-(£xH) 


Equation  B.  4 can  be  written  as 


f[n.(V-xH)e‘jl*,R  +r. 


Also,  from  Maxwell's  equation 


(n  x H)  jk  e"jkR  'RJ 


dS'  = o 
(B.  5) 


V'  x H = -jwtg  E 


in  a source  free  region. 
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Thus, 


$ (ft-  E)  (~ ju)< q)  e ^kR  ’ R • (n'x  H)  jk  e ^kR 


• (frx  H)  - jk  (fr.  E)J  e jkR  R 


dS*  = 0 


It  is,  therefore,  true  that 


ff  (nxH)-jk(n-  E)J  e-^  R ds'  = 0 (b.6) 


From  Eq,  B.l  and  B.6,  the  radial  component  of  E(x,  y,  z')  is  proved 
to  be  zero  at  far-zone. 

At  far-zone,  consequently. 


E(x, y, z)  = ff  [jco/i(nxH)0]tdS'  (5.56) 


APPENDIX  C 


PROOF  OF  EQ.  5. 92 


The  purpose  of  this  appendix  is  to  show  the  steps  involved 
in  arriving  at  Eq.  5.92  from  Eqs.  5.90  and  5.91.  The  same  steps  were 
also  taken  to  get  the  expression  (Eq.  5. 106). 

From  Eq.  5.90  and  Eq.  5.  91, 

aPn-l(cos  8>On-iW)  _ aPn+l  (cos  B)Pn4l<kR> 
(2n+  1)  • k • pn_  ^(ka)  (2n+  1)  • k • p ' (ka) 

(5.90) 

I1A(R)  = ‘max  sin  k,d  - R)  <5-91> 


e2r  ■ “K2  +fr> 

cR2 


where 


Therefore, 


d = a + h 


d U a I d .2 

I E„r  I. . (R) dR  = — f sin  k (d  - R)(k2  to  .(kR)dR 

R=a  2R  1A  (2n  + 1)  • k • pn_  ^(ka)  k dR2  n_1 


U a I 

max 

(2n+  !)•  k • pn+  (ka) 


d 

f sin  k(d-  R)(k2  + 
a 


-\2 

r 


)pn  + 1(kR)dR 


(C.  1) 
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Let 


I 


1 


sin  k(d-  R)(k2 


4'vi^® 


b 

= I sin  k(d-  RXk2 
a 


)Pn+1(kR)  dR 


(C.2) 


Changing  a variable  kR  to  z, 


kR  = z 


• . kdR  = dz 


f)2  a2 

— = k2  -A_ 

3R2  az2 


The  integral  Ij  becomes 


I 


1 


kd  -s 

/ sink(d  - z){k2  + k2  — — ) p (z)  -3? 
ka  az2  "-1  k 


kd 

kf  sink(d  - z)  p 
ka  n 


kd  2 

i(z)  dz  + k / sink(d  - z)  — 
ka  az2 


Pn_  1(z)  dz 


et 


kd 

hi  = k f sink(d  - z)  p (z)  dz 
ka  n~  1 


(C.  4) 


(C.3) 
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and 


I 


12 


tvu  2 

k / sink(d  - z) p (z)  dz 

ka  3z2  n"  1 


(C.  5) 


Then 


kd 

1^  = k f sink(d  - z)  ^(z)  dz  . 
ka 


Let 

Pn_  jCz)  = u sink(d-z)dz  = dv 

3 

3z“Pn-l(Z)  = du  cosk(d-z)  = v 


By  part  integration 


hi  - k 


r 9^n  1^ 

cosk(d  - z)  Pn_  x(z)  - J cosk(d  - z)  — — dz 


kd 


z = ka 


Integrating  the  second  term  oi  1^  by  part  again,  we  get 


kd  dp  (z) 

J cosk(a  - z)  — dz 

ka  dz 


p 
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Finally, 


E2R  llA^’  ^ “ 


U a I 


max 
2n+  1 


pn  - l<kd> 

A-  l'<ka) 


- cos  k(d  - a) 


pn  - l(ka> 

Pn-  !««) 


pn+  l<kd> 
pn  + l'<ta> 


Pn  , i(ka) 
+ cos  k(d  - a)  — * 1 

pn  + l’(ka> 


(5.92) 


APPENDIX  D 


NUMERICAL  EVALUATION  OF  Bn 
The  functions  shown  in  the  evaluation  of  the  coefficients  B 

n 

of  an  infinite  series  expressing  an  induced  surface  current  on  a hemi- 
spherical ground  plane  are  the  weighted  spherical  Kankel  functions  of 
the  second  kind  andtheir  derivatives. 

These  functions  Pn(x)  ci:n  be  expressed  in  the  form  of  a 
series  such  as 

p (x)  = xhj^(x)  = (j)n+1  e‘jX  Yj  (n  + -| , k)  — --j-  (D.  1) 

n n k=o  Z (2jxf 


where 


(n  + k)! 
k ! P(n  - k + 1) 


r(z  + 1)  = z!  z:  integers 

h^(x)  =<\^F  Hn+L(2)  (x)  and  j = ^ 

2 

Also,  from  the  recurrence  relations  of  Hankel  functions,  the 
following  relationship  between  the  spherical  Hankel  functions  and  their 
derivatives  can  be  obtained  {Ref.  18  ). 
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f h»w  = h {2)(x)--^-h(2)(x) 
dx  n n-1  x n 


(D.2) 


and 


P - 4 (*  hn<2)w)  ■ hf 1 'w + x -ar  hi2)«  <D-  3> 


dpn(x) 
dx 


From  Eqs.  D.  2 and  D.  3,  the  following  relationships  are 


derived 


P'0W 


d PQ(x) 
dx 


hf’fx)  - x hf(x) 


(D.4) 


P'n(x) 


d Pn(x) 
~dx 


x hn.j(2)(x)  - n h®(x) 


(D.5) 


p (x)  - x h (£,(x) 
*n  ' n 


(D.  6) 


Equations  D.  4,  D.  5,  and  D.  6 are  used  with  a table  for  spheri- 


cal Hankel  functions  to  evaluate  the  coefficients  B . 

n 


APP.  .NDIX  E 


EVALUATION  OF  R in  Eq.  5. 168 

rl 


From  the  power  series  representation  of  cosine  integral 
Ci(z),  Eq.  5.  165,  the  terms  associated  with  cosine  integral  of  a small 
argument  can  be  written  as 


lim  Ci(2k(a  + h)ej  * y + in  2k(a  + h)  + lim  in(e) 
e-o  ' ' c-o 


(E.l) 


lim  Ci(2kae) 
e-o 


lim  Ci  k(2  a + h)  e 
e-o 


y + ln(2ka)  + lim  fn(e ) 
6-0 


* y 


(E.2) 


+ in(k(2a  + h)V  lim  in(e)  (E.3) 

' / 6-0 


lim  Ci(kh6)  * y + in(kh)  + lim  in(e) 
e-o  e-o 


(E-4) 


Substituting  Eq.  E.  1 through  Eq.  E.  4 into  appropriate  places  in(T)through 

© 


R =60 
rl 


in  2 - lim  in  6 + y + lim  in(e)  + in(kh)  - cos  2k(a  + h)  in(2ka) 
e-o  e-o 


- cos  2k(a  + h)  in  2k(a  + h)  + cos  2k(a  + h)  in(2a  + h)k  - Ci(2kh) 


+ cos  2k(a  + h)  Ci(4ka)  + sin  2k(a  + h)  Si(4ka)  - (1  - ^^ka^  ^ s*n  ^ ^ 
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- cos  2k(a  + 
+ cos  2k(a 


sin  2k(a  + 


+ sin  2k(a  + h)  Si  ( 4k(a  + h) 


>i  ^4k(a  + h)  j 


(E.  5) 


It  should  be  noticed  that  the  coefficients  of  y and  in(e) 
associated  with  the  cosine  integrals  add  up  to  1.  This  can  be  shown 
easily  after  sflUTr  algebraic  manipulation  with  trigonometric  identities. 

After  cancelling  out  the  logarithmic  singularities  in  Eq.  E.  5 
and  collecting  terms,  the  final  form  of  Eq.  E.  5 can  be  shown  as  Eq.  5.  168. 
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